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ABSTRACT 

The LSPM-north catalog is a comprehensive list of 61,977 stars north of the J2000 celestial equator that 
have proper motions larger than 0.15 seconds of arc per year (local-background-stars frame). The catalog has 
been generated primarily as a result of our systematic search for high proper motion stars in the Digitized Sky 
Surveys using our SUPERBLINK software. At brighter magnitudes, the catalog incorporates stars and data 
from the Tycho-2 Catalogue and also, to a lesser extent, from the All-sky Compiled Catalogue of 2.5 million 
stars. The LSPM catalog expands considerably over the old Luyten (LHS, NLTT) catalogs, superseding them 
for northern declinations. Positions are given with an accuracy <100 mas at the 2000.0 epoch, and absolute 
proper motions are given with an accuracy of ?a8 mas yr"' . Corrections to the local-background-stars proper 
motions have been calculated, and absolute proper motions in the extragalactic frame are given. Whenever 
available, we also give optical Bj and Vj magnitudes (from TYCHO-2, ASCC-2.5), photographic Bj, Rp, In 
magnitudes (from USNO-Bl catalog) and infrared J, H, K, magnitudes (from 2MASS). We also provide an 
estimated V magnitude and y-7 color for nearly all catalog entries, useful for initial classification of the stars. 
The catalog is estimated to be over 99% complete at high Galactic latitudes {\b > 15|), and over 90% complete 
at low Galactic latitudes (\b > 15 1), down to a magnitude V = 19.0, and has a limiting magnitude V = 21.0. 
All the northern stars listed in the LHS and NLTT catalogs have been re-identified, and their positions, proper 
motions and magnitudes re-evaluated. The catalog also lists a large number of completely new objects, which 
promise to expand very significantly the census of red dwarfs, subdwarfs, and white dwarfs in the vicinity of 
the Sun. 

Subject headings: astrometry — surveys — stars; kinematics — solar neighborhood — stars; white dwarfs — 
stars; Population II 



1. INTRODUCTION 

The search for and identification of stars with large proper 
motions has traditionally been the dominant method for find- 
ing nearby stars. Stellar distances are ultimately determined 
with parallax measurements, but these generally require sub- 
stantial effort, and until the HIPPARCOS mission (1991), 
large systematic parallax surveys were impractical. One 
therefore had to rely on secondary diagnostics of proximity 
that could be more easily measured, such as large proper mo- 
tions. It is a historical fact that the vast majority of the stars 
within 25 parsecs of the Sun (the local volume of space known 
as the "Solar Neighborhood") have first been identified as 
high proper motion stars. These include some of our clos- 
est neighbors, such as Proxima Centauri (.Innes..l915,) . and 
Barnard's Star (Barnard 1916), which still is the star with the 
largest known proper motion (/i ~ 10.3" yr"'). 

Massive, large scale searches for high proper motion stars 
have been performed, and continually improved, throughout 
the 20th century, thanks to innovations in wide field astropho- 
tography . Early compilation s of known stars with large proper 
motions Jvan MaanenI 1 9 1 51) . were expanded by catalogs such 
as Max Wol f's Katalog von 1053 starker bewegten Fixternen 
(|Wolf'1919'), which he complemented over the years with sev- 
eral new additions, naming a total of 1567 stars after himself. 

' Based on data mining of the Digitized Sky Surveys, developed and op- 
erated by the Catalogs and Surveys Branch of the Space Telescope Science 
Institute, Baltimore, USA. 

- Developed with support from the National Science Foundation, as part 
of the NASA/NSF NStars program. 

^ American Museum of Natural History Kalbfleich Research Fellow. 



Frank E. Ross also contributed numerous discoveries, pub- 
lishing lists of new high proper motion stars over a period 
of 14 years (1925-1939), discovering a total of 1080 nearby 
objects (Ross 1939). These early investigations were made 
with the use of a visual blink comparator; photographic plates 
obtained at different epochs were blinked in succession, and 
examined by eye. The typical motions of the stars detected 
(which were at that time simply called "proper motion stars") 
were w 0.2- 1.0" yr-i. 

Over the years, the discovery of increasingly fainter stars 
having large proper motions pointed to the existence of 
a significant population of low luminosity stars. Some 
of the most extreme examples were discovered as faint, 
common proper m otion companions of brighter objects 
(|van Biesbroeck 1961). The very first free-floating brown 
dwarf was also discovered in a survey of high proper motion 
stars (Ruiz, Leaaett, & Allard 1997). Besides being an ex- 
tremely useful tool in the identification of nearby populations 
of low luminosity objects, proper motion surveys are also very 
sensitive to high-velocity stars. Because high-velocity (e.g. 
thick disk, halo) stars are detected out to larger distances than 
nearby disk stars in proper motion selected samples, they are 
significantly over-represented in catalogs of high proper mo- 
tion stars. Far from being a problem, this makes proper mo- 
tion catalogs highly useful tools for the study of Galactic stel- 
lar populations (including low-luminosity halo stars). This 
has motivated intensive searches for faint high proper motion 
stars over the whole sky. 

The largest deep surveys of high proper motion stars were 
carried out by the Lowell Observatory, and by Willem J. 
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Luyten at the University of Minnesota. Results from the Low- 
ell Proper Motion Survey were compiled and published in 
two lai-ge catalogs JGiclas. Burnham. & Thoma slT97llfT978i) 
listing a total of 1 1 ,749 stars with proper motion larger than 
0.2" yr"'. Luyten's work, on the other hand, was con- 
stantly updated over the decades, culminating in the publi- 
cation of two major catalogs: the New Luyten Catalogue of 
stars with proper motions larger than two tenths of an arc- 
second (Luvten 1979b) known in short as the NLTT, listing 
58,845 objects, and A catalogue of stars with proper motions 
exceeding 0".5 annually (Luvten 1979aJ known as the LHS, 
and which is essentially a subset of the NLTT listing 4470 of 
the fastest-moving stars. 

Most of Luyten's success stems from his use of the Na- 
tional Geographic Palomar Sky Survey (POSS-I), completed 
in the 1950s. Luyten obtained second epoch images in the 
late 50s and 60s using the same instrument and setup (at the 
Palomar Schmidt telescope). Luyten also developed a laser- 
scanning microdensitometer machine to process the north- 
ern sky images, considerably improving over previous eye- 
blinking methods. The depth of the Palomar plates allowed 
him to probe deeper than anyone before, mapping the turnover 
in the lo cal luminosity function at the bottom of the main se- 
quence jLuyteniri968h : he was however limited by the depth 
of his second-epoch plates, which were about 1 magnitude 
shallower (19th magnitude Hmit) that the POSS-I plates (20th 
magnitude limit). Until today, the NLTT catalog remains the 
largest and most complete list of high proper motion objects, 
at least for declinations north of -32.5° (the southern limit of 
POSS-I). 

Because the LHS and NLTT catalogs contain large num- 
bers of astrophysically significant objects, they have been 
used as a source of targets in many follow-up programs, 
with the faster LHS stars naturally taking precedence. Pho- 
tometric studies have included the search for nearby dwarfs 

tfeis 1984, 1996), multiband studies of halo stars ( Rv^ 
89) and of low-luminosity dwarfs and subdwarfs ( BesseUl 
91.) . Spectroscopic follow-up surveys have resulted in 
the identification of ne w nearby stars (Gizis & Reidlll99'7h . 
cool halo subdwarfs (Ruiz & Aneuita 1993), and white 
dw arfs (Hintzen. P. 1986; Bergeron. Ruiz. & Leggett 199^ 
IVennes & Kawkall2003l) . Despite numerous studies and ob- 
serving programs, hundreds of LHS stars and the majority of 
NLTT stars are still lacking formal spectral classification. The 
NLTT catalog in particular remains a goldmine of astrophysi- 
cally interesting but uncharacterized targets. Perhaps the most 
intensive follow-up study to date is a recent program devoted 
to the identification of nearby stars missing from the census 
of objects within 25 par secs of the Sun (Reid & Cru zll2002t 
jReid^ Kilkenny. & Cruzll2002HCruz &'Reidiii2002£ iReid et al\ 

Modern astrometric techniques require much more accu- 
rate positions and proper motions than initially recorded 
by Luyten. One example i s the identification of future 
possible microlensing events iSalim & Gouldll2000t) . Mo- 
tivated by these requirements, proper motions and posi- 
tions of the LH S stars have been recently recalculated by 
iBakos. Sahu. & Nemeth (2002). who systematically searched 
for all the stars in the Digitized Sky Surveys. A revision 
of the NLTT positions and proper motions was also initi- 
ated by Gould &Salim (2003) and Salim & Gould (2003), 
using the USNO-A2 catalog (Monetefa/. 1998) as a first 
epoch, and the 2MASS Second incremental Release as a sec- 
ond epoch. These revisions revealed that a significant fraction 



of LHS/NLTT entries contained large positional errors, up to 
several minutes of arc in some cases. Indeed, a comparison 
of NLTT and LHS positions, for stars appearing in both cata- 
logs, has revealed the existence of ty pographical errors in both 
catalogs (Lepine. Shara. & Rich 2002). The uncertainties in 
Luyten's positions explain in part the difficulty in carrying out 
follow-up observations of his objects, and raises the possibil- 
ity that background stars have been mistaken for NLTT stars, 
resulting in erroneous spectral classifications. This motivates 
a complete re-evaluation of the positions of all NLTT objects. 

The problem of the completeness of the NLTT has been 
much debated (see Pokorny, Jones, & Hambly 2003, and 
references therein). The reason is that statistical studies of 
the Luyten stars (LHS stars in particular) have been used 
in estimates of the local density of low-mass and degener- 
ate stars. In particular, the NLTT has been used to estimate 
the local density and luminosit y function of white dwarfs 
(iLieb ert. Dahn. & Mo net' '1988^, 'Leggett. Ruiz. & Be rgeronI 
1998) and the density and luminosity function of the local 
halo population ("Dawson 1986; Lee 1991; Gould 2003a). The 
accuracy of those statistical studies, however, is entirely de- 
pendent on the completeness of the underlying sample, or at 
least on their estimated completeness. Hence the importance 
of being able to estimate the completeness (as a function of 
position, magnitude, proper motion, ...) of the NLTT catalog. 

As initially noted by Dawson ( 1986), the NLTT and LHS 
catalogs are notably incomplete in two distinct areas: (1) 
south of -32.5° in declination, and (2) in a band within ±10° 
of the Galactic plane. However, outside of these specific ar- 
eas, in an area referred to as the Completed Palomar Region 
(CPR), the catalogs do appear to be significantly complete 
(> 80%) at least down to 19th magnitude. This is significant 
because the magnitude distribution of NLTT stars reflects the 
turnover at the faint end of the luminosity function: the num- 
ber density of NLTT stars peaks at /? « 15, well before the 
limiting magnitude of the POSS-I plates on which they are 
based. Because they probe the luminosity function turnover, 
and because they contain significant samples of both disk and 
halo stars, the LHS and NLTT catalogs are major tools in 
the determination of the local luminosity function of both the 
Galactic disk and Galactic halo. In any case, the accuracy of 
the luminosity function and of the stellar density is dependent 
on a proper evaluation of the completeness of these proper 
motion selected samples (Gould 2003a). 

The completeness of a proper motion catalog can be deter- 
mined eith er internally , or ex ternally. One internal test was 
devised by Flvnn et al\ ll200li) . and uses the fact that both the 
magnitude and proper motion are a function of distance. For 
example, if one takes subsamples with a lower proper motion 
limits fix and /i2 = 1.259/ii, sample 1 should contain objects 
that are on average more distant by a distance modulus of 0.5. 
The completeness at a magnitude V\ in sample 1 can thus be 
determined by comparing the relative number of stars in bins 
V\ and V2 = Vi -0.5 in sample 2. One can then work iteratively, 
starting from a bin V2 assumed to be complete. Applying the 
technique to the NLTT, Flynn et al. (2001) estimated that the 
completeness in the CPR (for stars with ji > 0.2" yr"') falls 
linearly from 100% at V = 13 to 60% at V = 18, and then 
breaks dow n to 0% at V = 20. T he method was, however, put 
in doubt by [Monet et al.\ ( l2000h who noted that at high galac- 
tic latitudes the density of stars decreases with distance, and 
sample 2 is thus not equivalent to the more distant sample 1 ; 
the method of Flynn et al. (2001J would thus tend to underes- 
timate the completeness. iGouldl (l2003a^ discusses this prob- 
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lem (see his appendix), a nd conc ludes that there is little evi- 
dence to suggest the Flv nn et al\ (llOOl) result to be in error 
The main problem is a lack of reliable external completeness 
estimates. 

External completeness tests are based on a direct compari- 
son of the proper motion catalog (the NLTT) to a deeper or 
more sensitive proper motion survey conducted over a se- 
lected area. The main caveat of that method is that the re- 
sults then become dependent on the completeness of the new 
survey itself, which has to be estimated by other means. Of 
course, if the new survey is conducted over a sufficiently large 
area, it can simply supersede the former catalog. 

From a small sample of only 100 square degrees from the 
APM Proper Motion Project, Evans ( 1992) concluded that the 
NLTT has an incompleteness of w 16%, most of it probably 
due to random measurement errors at the lower proper mo- 
tion limit of 11 = 0.18" yr"'. This compl eteness level is larger 
than suggested by the internal test of Fl vnn et al\ (1200 h) . but 
the small area and limited number of objects puts this exter- 
nal test in doubt. A much larger survey of 1,378 square de- 
grees in the northern sky (Monet et al. 2000) was performed 
using pairs of pl ates from the Second Palomar Sky Survey 
dReid et a/.lll991ii) . Because of the short timespan between 
pairs of plates, this survey was limited to the detection of 
stars with very large proper motions (0.4" yr"'). Neverthe- 
less, fifteen stars were found that had been missed by Luyten, 
of which 6 have proper motions ji > 0.5" yr~', suggesting a 
completeness > 90% for the LHS catalog. (N.B. while the 
Monet et al. paper cite s 17 "n e w" obj ects, two of them were 
subsequently found by Gould (2003a) to be the NLTT cata- 
log stars 58785 and 52890.) Unfortunately, their survey did 
not cover the 0.18" yr"' < /x < 0.4" yr"' range, where most 
of the NLTT stars are found, leaving open the question of the 
NLTT completeness. 

Many more attempts have been made at conducting 
more sensitive proper motion surveys in the south. Be- 
cause the NLTT is well known to be much less complete 
south of -32.5°, the potential pay-off is much higher. In 
the Calan-ESO Proper Motion Survey - famous for its 
discovery of Kelu-1, the first free floating brown dwarf 
(Ikuiz. Leggett. & Allard 1997) - 14 pairs of ESO Schmidt 
plates were used , covering an area of 350 square degrees 
dRuiz et al\ Fourteen new stars with ji > 0.5" yr ' 

were found, suggesting a completeness of only sa 60% for 
the LHS in the south. But since several of their areas 
are south of -32.5°, this overestimates the incompleteness 
of the NLTT in the CPR. Th e larger survey co nducted by 
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IWroblewski & Torres (1989, 199 J 
continued by Wroblewski & Cost^ (119991 
ers a total of 3,275 square degrees in 131 scattered areas. 
This survey is performed by direct visual inspection of pho- 
tographic plates, using a Zeiss-Jena plate comparator The 
survey now has 2,495 cataloged objects, all of which are new 
(i.e., not in the NLTT). Within the limits of the CPR, they 
typically discover «50 new stars with /i > 0.2" yr"' in every 
100 square degree area. When compared to the NLTT density 
of roughly 90 stars per 100 square degrees in the same areas, 
this yields an overall estimated completeness of « 65% for 
the NLTT, but a more detailed analysis would be warranted. 

Other large southern surveys are built from lists of ob- 
jects generated from machine scans of photographic Schmidt 
plates. A survey of 2,000 square degrees near the south polar 
galactic cap was made with lists of objec ts from UK Schmid t 
plates, scanned with the APM machine dScholz e?a/.ll200fll) . 



This survey is ongoing, and may eventually cover much larger 
areas of the southern sky. The Liverpool-Edinburgh h igh 
proper motion survey ("Pokor nv. Jones. & Hamblvl l2003h is 
based on lists of objects from ESO Schmidt and UK Schmidt 
plates, scanned with the SuperCOSMOS machine, and has a 
lower proper motion limit of 0.18" yr"'; the survey is cur- 
rently limited to a moderately large area of the southern sky, 
covering « 3,000 square degrees around the south Galactic 
cap. More recently, Hamblvef a/. (2004) have used Super- 
COSMOS data too, to search w 3,000 square degrees south 
of -57.5° for stars with proper motion /i > 0.4" yr"' (although 
they only list those they found having > 1 .0" yr"'). What is 
clear from these southern sky surveys is that the LHS/NLTT 
catalogs are definitely incomplete south of -32.5°. The num- 
ber of new high proper motion stars discovered is a sizable 
fraction of all the high proper motion stars identified. The 
consensus is that in the southern sky, the NLTT is no more 
than w 70% complete down to R=19, and possibly much less 
complete in some areas. 

More massive astrometric catalogs based on lists of objects 
from scanned Schmidt plates are in preparation, including the 
second Guide Star Catalog (GSC-11). By combining lists of 
objects from multiple epochs, these will attempt to provide 
proper motions for al l stars detected. T he already available 
USNO-Bl.O catalog (Monet ef fl/."2003) is a first attempt at 
such a deep, all-sky, astrometric catalog with proper motions. 
For the northern sky, the plate material includes the Oschin 
Schmidt plates from the first and second epoch Palomar Sky 
Surveys (POSS-1, POSS-11), providing a large temporal base- 
line and deep (R=20) limiting magnitude; in the southern sky 
scans of the ESO schmidt and UK schmidt were also used. 
Unfortunately, the huge number of objects involved (1 billion 
sources), and the large number of false detections (plate de- 
fects, mismatches from multiple epoch detections) makes the 
identification of high proper motion stars very difficult. An 
analysis of the USNO-Bl.O catalog by Gould (2003b) shows 
a huge contamination rate (200 to 1) in bo gus high proper 
motion objects. Recent efforts (^Gould & Kollmeied 120041: 
iMunn e t al. 2004) have however shown that it is possible to 
eliminate most of the false detections by cross-correlating the 
USNO-Bl.O catalog with other large area surveys like the 
Sloan Digital Sky Survey (SDSS). Unfortunately, the com- 
pleteness of the USNO-Bl.O for high proper motion stars 
(based on the recovery of NLTT stars) does not exceed 90% 
over the whole sky, and drops to 70% at low galactic latitudes. 

Very high completeness levels have however been achieved 
by us, in a novel approach to surveys of hig h proper mo- 
tion stars ( Lepine. Shara. & RichI 120021 120031) . Instead of 
working with lists of objects identified from plate scans, 
we work directly with the pixel data, using an image sub- 
traction algorithm. With the help of a specialized soft- 
ware (SUPERBLINK), we have performed a very success- 
ful proper motion survey based on a massive re-analysis 
of all image scans of the Oschin Schmidt plates (POSS-I 
and POSS-II) made at STScI for the Digitized Sky Surveys 
(DSS). The image subtraction method is significantly more 
efficient in densely populated fields, where examination with 
the blink comparator is difficult, and identification of point 
stellar sources by scanning machines considerably less effi- 
cient because of crowding. To massively apply plate subtrac- 
tion methods to the DSS data, we have developed the SU- 
PERBLINK software, which works like an automated blink 
comparator. Our initial survey was a search for stars with very 
large proper motions /i > 0.5" yr"' over the whole 20,000 
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square degrees of northern sky. We showed our method to be 
extremely successful, by recovering essentially all LHS stars, 
and with the discovery of 198 new stars with /i > 0.5" yr"'. 
This yielded the definitive completeness measure of the LHS 
catalog in the northern sky (1662/1866= 89.1%) and con- 
firmed the results of Monet et al, ( 2Q0 0) about the northern 
sky completeness of the LHS. However, it left open the ques- 
tion of the completeness of the NLTT. Our new proper mo- 
tions catalog directly addresses this problem since it is an ex- 
tension of our survey to smaller proper motions (ji > 0.15" 
yr"'), with a completeness level exceeding 99% down to 
Rp = 19. The executive summary is: in the northern hemi- 
sphere, we find and report 40,843 stars with /i > 0.18" yr~^ , of 
which 28,486 are NLTT objects; the LSPM catalog also lists 
an additional 21,133 stars with 0.15" yr"' < fi < 0.18" yr"', 
only 2,875 of which are in the NLTT catalog. Our analysis 
demonstrates that in the CPR, the completeness of the NLTT 
is 85% at V = 18, and breaks down at V = 19, while at low 
galactic latitudes (\b\ < 10) its completeness falls from 90% 
at V = 15 to only 30% atV =11. This suggests that the internal 
test of Flvnn et al. (2001) underestimates the completeness of 
the NLTT in the CPR. 

Note that other major proper motion surveys, which have 
been extremely successful in determining highly accurate 
proper motions of selected stars, are not very helpful in in- 
creasing the completeness of our proper motion catalogs. 
These include the astrometric survey conducted with the 
HIPPARCOS satellite, w hose data ar e now compiled in the 
TYCHO-2 catalog (Hog ef aZ.ll2000h . Only stars brighter 
than V=9 were observed systematically, and an input cat- 
alog was used for stars down to the limiting magnitude 
(V=13). Bright NLTT stars were included in the survey, 
but very few new high proper motion stars were discov- 
ered above Luyten's cutoff. Also limited as tools for finding 
new high proper motions stars are the Lick Northern Proper 
Motion Program ( Klemola. Jones. & Hanson 1987) and its 
southern extension, th e Yale/San Juan Proper Motion Survey 
dvan Altena & L6pe5iJ99 1). Both programs aim at very pre- 
cise astrometric measurements of selected stars, and largely 
rely on an input catalog, although a subset of stars were 
picked at random. Another highly accurate astrometric sur- 
vey is the US N aval Observatory CCD Astrograph program 
dZacharias et ah .2000). from which an all-sky astrometri c 
catalog (UCAC) is being assembled dZacharias et aZ.ll2()03l) . 
While extremely promising as an expansion to the TYCHO-2 
catalog, the UCAC will have a limiting magnitude of R=16, 
and it is unclear how sensitive the survey will be for stars with 
very large proper motions. At this time, it appears that our im- 
age subtraction method holds the best promise for generating 
an all-sky replacement to the LHS/NLTT catalogs. 

As part of the NASA/NSF NStars initiative, we have been 
expanding our DSS-based survey, aiming at the systematic 
detection and verification of all stars in the northern sky with 
proper motions larger than 0.05" yr"'. Our goals are to 
achieve optimal detection rates, with completeness exceeding 
99% down to 7? = 19.0 over most of the sky, with minimal con- 
tamination from false detections. This paper presents our first 
major data release: a catalog of all known stars in the northern 
sky with Proper Motions larger than 0. 15 " yr"' . This catalog 
which we refer to as the LSPM catalog, includes improved 
astrometry and photometry for more than 3 1 ,000 high proper 
motion stars previously listed in the LHS and NLTT catalogs. 
The LSPM catalog also incorporates bright high proper mo- 
tion stars from the TYCHO-2 catalog. Finally, the LSPM cat- 



alog contains over 28,000 newly discovered high proper mo- 
tion stars. 

The LSPM catalog represents a major improvement over 
the NLTT catalog. Not only is it much more complete, but the 
positions and proper motions are also much more accurate. 
In effect, the LSPM supersedes the NLTT for the sky north 
of the celestial equator, and should now be used in its place 
for all applications. This paper provides information that is 
essential in understanding how the LSPM catalog was built, 
and what are its strengths and limitations. A description of 
the SUPERBLINK code, used to find LSPM candidates in the 
DSS, follows in §2. The detailed procedure for the inclusion 
of a star in the LSPM catalog is detailed in §3. The sources 
used for the photometry are presented in §4, while the catalog 
astrometric accuracy is discussed in §5. The format of the 
catalog is explained in §6. The completeness of the LSPM is 
discussed in §7. A preliminary analysis of the stellar contents 
of the catalog is given in §8. Plans for future expansion and 
improvement of the catalog are summarized in the conclusion 
(§9). 

2. A NEW SURVEY FOR HIGH PROPER MOTION STARS 

2.1. The SUPERBLINK software 

SUPERBLINK is an automated bl ink comparator devel- 
oped by SL, and first described in 'Lepine, S hara. & Ricl] 
(2002). Given two different images of the same patch of sky 
on input (up to 2kx2k pixels in size), SUPERBLINK auto- 
matically identifies any object that has moved between the 
two epochs, such as a star with a large proper motion. On out- 
put, the software generates a list of possible moving objects 
in the field, with their positions, proper motions, image mag- 
nitude, and a probability index that estimates the likelihood 
of the object being real. The software also generates iden- 
tification charts (151 x 151 pixels in size) centered on each 
object. These charts are dual-epoch, and can be blinked on 
the computer screen for easy examination of the moving ob- 
ject. The two core elements of SUPERBLINK are an image 
superposition and subtraction algorithm (SUPER), and a shift- 
and-match search algorithm (BLINK). 

The current version of the code has been optimized for use 
with Digitized Sky Survey images (first epoch DSS and sec- 
ond epoch XDSS) to look for stars with large proper motions. 
Pairs of images are provided to the code on input, each pair 
consisting of one 17' x 17' field extracted from the DSS at a 
specified position, and a second 17' x 17' extracted from the 
XDSS and centered on the same position. In the northern sky, 
the DSS image invariably consists of a POSS-I scan, with a 
typical resolution of 1.7" per pixel, while the XDSS image 
is a POSS-II scan with a resolution of 1.01" per pixel. The 
POSS-II image is generally of higher quality than the POSS- 
I image: the background noise is lower, the image reaches 
about one magnitude deeper, and the astronomical resolution 
(seeing) is better. 

The SUPER procedure performs image transformations to 
make the two images in the pair look as similar as possible. 
The procedure uses the first image (DSS) as a template, and 
attempts to modify/degrade the second image (XDSS) in such 
a way that it can be subtracted from the first with the smallest 
residuals possible. The SUPER procedure follows a series of 
steps described below. 

Rescaling. The two images, on input, can have different res- 
olutions; the SUPER procedure resets the two images to the 
same angular scale. Typically, the higher resolution second 
epoch image (XDSS) is remapped onto a grid that matches 
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the resolution of the lower resolution first epoch (DSS). 

Rectification. Images are rectified so that their background 
levels (sky) are uniform and set to a value of 1. The code 
uses a procedure that marks each pixel as either "sky" or 
"object" (based on the statistics of intensity values). A two- 
dimensional linear fit is then performed on the "sky" pixels. 
The image is then divided by this fit, setting background levels 
to unity. While the background level is never strictly linear on 
a photographic plate (edge effects are important), it is a good 
approximation on the scale of the images provided on input 
(17' X 17'), which are much smaller than the typical size of 
the POSS plates (384' x 384'). 

Normalization. The total flux from all "object" pixels is de- 
termined for each image. The second image is then normal- 
ized so that the total flux (above background) from "object" 
pixels is equal to that in the first image. Note that this normal- 
ization might be inaccurate if there are bright objects showing 
up in only one of the two images. This does happen, particu- 
larly if the two images are not exactly aligned initially. It may 
then happen that e.g. a bright star near the edge of one image 
does not show up in the other, and vice versa. The renormal- 
ization procedure (see below) will generally correct for any 
normalization error. 

Shift: and rotate. The second image is shifted vertically and 
horizontally (AX, AY), and rotated (A0), before being sub- 
tracted from the first image. The procedure is repeated recur- 
sively, first using small incremental values of AX, AY, and 
A6' until a good match is found, i.e. until the residuals are 
significantly smaller than the total flux in each image. More 
precise values for the shift and rotation are then determined 
using a multidimensional downhill simplex minimization rou- 
tine, which identifies a strong minimum in the residuals in 
[AX, AY, A6'] space. Note that this procedure accurately su- 
perposes the two images using the assumption that most stars 
in the field are "fixed". Any systematic motion in the back- 
ground "fixed" stars will be eliminated. This means that all 
proper motions calculated by SUPERBLINK will be proper 
motions relative to the background of 'fixed" stars, and not 
absolute proper motions. 

Renormalization. The second image is normahzed again, 
as described above, but this time using only "object" pixels 
that are common to both images in the pair. These can now 
be easily determined since we know from the preceding shift 
and rotate procedure which part of the field is common to both 
images. 

Convolution. The second image is then degraded so that its 
PSF matches that of the first one. A convolution profile of 
variable width is applied to the second image, which is then 
subtracted from the first image. The width of the profile is in- 
creased until a minimum value in the residuals is found. The 
shape of the convolution profile has been determined by trial 
and error. Several different shapes have been considered; a 
simple profile generated by a sum of two gaussians of differ- 
ent widths was found to yield the best results. This same gen- 
eral profile was applied to all our fields. After this final pro- 
cedure, the first and second images generally look extremely 
similar. In the best of cases, it is very difficult to tell the two 
images apart just by looking at them. The only obvious differ- 
ences are variations in the noise patterns, or the presence of a 
variable star, an asteroid track, or a high proper motion star. 

The BLINK procedure starts with one pair of images that 
have first been superposed with SUPER, and proceeds to iden- 
tify any object that has moved between the first and second 
epoch. Stars with very large proper motions essentially appear 



as pairs of objects, one at each epoch, that do not cancel out 
after image subtraction; these are fairly easy to find. On the 
other hand, stars whose total motion between the two epochs 
is less than their apparent sizes on the POSS plates show a 
more complex pattern in the residuals, having been partially 
canceled out. On scans of POSS-I plates, typical sizes (full 
width at half maximum) of stars range from « 3 pixels (5.1") 
for unsaturated stars (R/ > 15), to w 15 pixels (25.5") for 
the brightest, saturated objects detectable by SUPERBLINK 
(Ry K, 10). The minimum motion of stars in our catalog is 
about 6" between flie POSS-I and POSS-H plates (0.15" yr"' 
in 40 yr). This means that while the fainter stars are always 
well-separated after plate subtraction, the images of many of 
the brighter proper motion stars will overlap, and will par- 
tially cancel out after subtraction. The following procedures 
in BLINK allow for a correct treatment of all moving objects, 
whether or not they partially cancel out on plate subtraction. 

Subtraction and cataloging of residuals. The two images 
processed by SUPER are subtracted from each other. Any 
object that has moved significantly between the two epochs 
induces a large, local maximum/minimum on the residual 
image. All the minima/maxima are mapped, cataloged, and 
matched to their source on the first or second image. 

Removal of the candidate moving objects from the first 
epoch image. Each object on the first epoch image that is asso- 
ciated with a large residual is removed from that image, with 
all pixels set to the sky level. A new residual image is then 
calculated by subtracting the second image from the first one 
(which now lacks the profile of the candidate moving star). 
Because the moving object has now been removed from the 
first image, its second epoch profile now shows up in its en- 
tirety in the residuals. This allows for an easy identification of 
the slow-moving objects, which would otherwise be canceling 
themselves out partially on the residual image. 

Search, and match the moving star on the second epoch. 
The profile of the object that was removed from the first image 
is then recursively shifted in X and Y, and added to the resid- 
ual image. If the star is indeed a moving object, there will be 
a near exact replica of it on the residual image. Because the 
second epoch has been subtracted from the first, the rephca 
will be a negative source of similar flux. Hence, the object 
will cancel out its second epoch replica when it is shifted by 
the AX and AY which corresponds to its motion on the plate 
between the first and second epoch images. 

Calculate likelihood. For each candidate moving object se- 
lected on the first epoch image, the code identifies the best 
possible match for that star on the second epoch image, within 
a radius of 1.5'. The software distinguishes between actual 
moving objects, and accidental matches of unrelated features 
based on the quality of the match. A probability index is cal- 
culated for each candidate moving object which is function 
of: (1) the difference in the object magnitude between the 
first and second epoch, and (2) the difference in the magni- 
tude density of the object between the first and second epoch. 
A detection thus has a high likelihood if the object has the 
same magnitude on both images, provided that the type of 
object (compact/extended) is also the same. Criterion (2) es- 
sentially prevents stars from being matched with galaxies, and 
vice-versa. 

Repeat, for candidate moving objects from the second 
epoch image. The procedure is repeated, but this time only 
the second epoch counterparts are considered. The object is 
subtracted from the second epoch image, residuals are recal- 
culated, and the object is shifted until it cancels out its first 
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epoch profile, now showing in its entirety in the residuals. 
Note that this means that most high proper motion stars are 
identified twice by the code, once from their first epoch loca- 
tion, and once from their second epoch location. This redun- 
dancy is necessary for the identification of stars whose profile 
is superposed on the profiles of other stars at either epoch (es- 
pecially in fields with significant crowding). In effect, this 
increases the chance of detection for blended stars; the star 
will be detected even if it blended with another source in ei- 
ther of the two epochs. This is especially useful for faint stars 
moving in the vicinity of brighter objects or in crowded fields. 

Once the image has been completely analyzed and 
searched, the code generates a list of all candidate moving 
objects along with their positions, relative proper motions (in 
pixels per year), integrated plate magnitudes, and likelihood 
index. The code uses the plate solutions and epochs (found in 
the image headers) to determine the local scale and orienta- 
tion of the plates, and calculate the magnitude (in seconds of 
arc per year) and direction of the proper motion vector on the 
sky. The software also extracts 151 x 151 pixels^ finder charts 
from the superposed images. These charts are extremely use- 
ful, as they are used to subsequently verify each and every 
detection by eye, on the computer screen. 

2.2. Application to the Digitized Sky Surveys 

The first epoch of the Digitized Sky Surveys (DSS) in the 
northern sky consist of scans of photographic plates from 
the National Geographic Palomar Observatory Sky Survey 
(POSS-I), obtained circa 1950. The scans were performed 
with the GAMMA machine by the Catalogs and Surveys 
Branch at the Space Telescope Science Institute. Only the 
red plates (xxlOSaE emulsion + plexi) have been scanned. 
The second epoch of the Digitized Sky Surveys (XDSS) con- 
sists, for the northern sk y, of scans of the Second Palomar 
Observatory Sky Survey (Reid et al. 1991). Images from the 
POSS-II include scans of plates from all three photographic 
bands of the survey: the blue (IllaJ emulsion with GG385 
filter), red (IIIaF emulsion with RG610 filter), and near in- 
frared (IVN emulsion with RG9 filter). The data from both 
the DSS and XDSS are publicly available from a variety of 
on-line databases 

We divided the northern sky into 615,800 areas distributed 
on a grid with a separation of 12' in DEC and a mean sepa- 
ration of 10' in RA. At every grid point, we extracted from 
the DSS (our first epoch) and XDSS (our second epoch - red 
band only) pairs of images each 17' x 17' on a side. We de- 
liberately extracted images that are much larger than the grid 
point separation, thus allowing for a significant overlap be- 
tween neighboring image pairs. 

We allowed for a large overlap between neighboring sub- 
fields in part because of the required rotation of one of the im- 
ages in the superposition process. Square subfields extracted 
from the DSS are generally not oriented with the Y axis point- 
ing toward the celestial north pole; rather they follow the lo- 
cal XY coordinates of the scanned POSS plates. As a result, 
pairs of images extracted from the DSS and XDSS are gen- 
erally not aligned, and the XDSS image is rotated (by up to 
30 degrees at high latitudes) by the SUPER procedure. Areas 
near the corners of the square subfields are thus cut out. We 
therefore allow for a band 1' wide running along the edge of 
each subfield so that no gaps in sky coverage occur 

A large overlap is also required for completeness because 

Including |http://archive.stsci.edu/| 



a high proper motion star, to be detected by SUPERBLINK, 
must be present in a given subfield at each of the two epochs. 
A star that has moved from one subfield to another would not 
be detected as a moving object but rather as two distinct "vari- 
able" stars. A star with a proper motion p < 2.0" yr"' can 
move up to 1.5' between the two epochs of the POSS-I and 
POSS-II. This is why we also allocated an additional band 
1.5' wide running along the edge of each field, to help in the 
detection of stars with very large proper motions. 

In summary, the different angular scales, scanning resolu- 
tion, non-alignment of subfields, different pixellation grids 
and offsets between scans, and different image quality and 
limiting magnitude are all accounted for and corrected by SU- 
PERBLINK. 

All of our DSS scans were extracted from The Digitized Sky 
Surveys series of CD-ROMs, published by the Space Tele- 
scope Science Institute (STScI). All the XDSS scans were 
downloaded off the Internet directly from the STScI archive 
(where they are stored on a CD-ROM jukebox), with kind 
permission of the STScI Catalogs and Surveys Branch. All 
subfields were processed as they were downloaded. Complete 
uploading/downloading and analysis of all 615,800 subfields 
was performed over a period of 11 months, from May 2001 
through March 2002. Computations were performed on a dual 
Pentium-Ill processor machine running Linux. Scripts were 
used to automate the procedure, and the downloading and pro- 
cessing of the whole northern sky with SUPERBLINK was 
completed with minimal user interaction. Most of our human 
effort went into the quality control phase, described in the next 
section. 

2.3. Visual confirmation of candidates 

False detections are inevitable when one is looking for high 
proper motion stars on photographic plates. The POSS plates 
are filled with plate defects of different sorts, such as grains 
and bubbles in the emulsion, dust specks, and scratches. The 
plates also contains transient images left by solar system bod- 
ies (asteroid tracks), and the occasional meteor trail, nar- 
row artificial satellite track (POSS-II only), or wide airplane 
track. A combination of plate defects and/or space junk may 
conspire to create the illusion, on a DSS/XDSS pair, of an 
object moving at a rate within our detection limits (0.15" 
yr-> <p< 2.00" yr"'). Of course, SUPERBLINK automati- 
cally eliminates most such bogus detections with the require- 
ment that candidate high proper motion stars must have com- 
parable fluxes and flux densities on both plates. However, 
it is not uncommon to see plate defects of the same magni- 
tude within seconds of arc, or minutes of arc of each other on 
two different epochs, mimicking the behavior of a high proper 
motion star. This is especially true for faint features near the 
detection limit of the plate, which tend to be very numerous. 

Another major source of false detections is the long diffrac- 
tion spikes associated with the brighter stars. Because fields 
from the DSS and XDSS often do not have the same orien- 
tation on the sky, the position angles of the diffraction spikes 
change from the first to the second epoch. Once the two im- 
ages are superposed and subtracted out in SUPERBLINK, 
diffraction spikes systematically show up as intense resid- 
ual features. When the superposed images are blinked on 
the computer screen, the spikes display a remarkable rotat- 
ing motion between the two epochs. This motion is of course 
recorded by SUPERBLINK which systematically lists mov- 
ing spikes as possible proper motion objects. One solution 
that was considered at first was to reject any detection of a 



Catalog of Northern Stars With /j, > 0.15" yr"i 



moving object within a certain distance of a bright star How- 
ever, after the detection of several faint high proper motion 
stars in the vicinity of bright diffraction spikes, we decided to 
investigate them all, to maximize the detection rate of genuine 
high proper motion stars. 

The most direct and reliable way to eliminate false detec- 
tions is by visually inspecting each and every candidate high 
proper motion star, using a blink comparator A trained eye 
can easily distinguish real stars from plate defects, for objects 
down to a magnitude of r » 19. The second epoch of the Dig- 
itized Sky Surveys (the XDSS) also contains images in the Bj 
and In band, which can be used as a third epoch for confirma- 
tion of ambiguous objects. 

Blinking each and every object identified by SU- 
PERBLINK is a daunting task. However, the task is actu- 
ally made easy (if only time consuming) thanks to the conve- 
nience of the finding charts generated by SUPERBLINK. The 
SUPERBLINK charts are more than just pairs of DSS/XDSS 
scans. While the first epoch of the chart is essentially the 
DSS image centered on the candidate moving object, the sec- 
ond epoch of the chart (as explained in §2.1 above) is an 
XDSS image that has been processed and modified by SU- 
PERBLINK to match the appearance and quality of the DSS 
image. Using simple software (designed by SL), it is possible 
to blink sequentially large numbers of finder charts, accepting 
and rejecting stars with a single keystroke, and automatically 
updating the list of confirmed high proper motion stars. With 
a little training, it is possible to quickly sift through hundreds 
of candidates, at a rate of about 1 star per second. False de- 
tections typically outnumber real objects by a factor of 3 to 4. 
The visual confirmation of w60,000 high proper motion stars 
carried out for this catalog thus represents a total of about 75 
hours of intensive human inspection. 

An interesting benefit of the 60,000 individual visual in- 
spections was the identification of close proper motion pairs. 
The SUPERBLINK software does not discriminate between 
point sources and extended objects, and in the course of the 
survey, several extended objects were found to be moving. 
Many of these turned out to be double stars with small sepa- 
rations (« 1-5"). Close pairs, on the POSS plates, produce 
images that are sometimes elliptical, if the two stars are of 
equal magnitude, and sometimes "pear-shaped", if the stars 
have different magnitudes. While these were all identified as 
single moving objects by SUPERBLINK, they were flagged 
as probable multiple systems during visual inspection. All 
those that could be confirmed were then listed as distinct ob- 
jects in the catalog (see §3.2 below). 

3. BUILDING THE LSPM CATALOG 

3.L Stars identified with SUPERBLINK in the Digitized Sky 

Surveys 

Using SUPERBLINK, we have successfully analyzed 
DSS/XDSS fields covering 99.23% of the northern sky 
(20,460 square degrees). Areas that were not analyzed in- 
clude a small patch of sky north of 87 degree in declination 
(» 30 square degrees), which we avoided because of prob- 
lems associated with the very large rotation angles required 
in the superposition of the first and second epoch image. SU- 
PERBLINK also failed to analyze some 4,766 scattered sub- 
fields, covering a total area of 295 square degrees; these were 
rejected after SUPERBLINK was unable to superpose the two 
images because of the presence of a very bright, saturated ob- 
ject in the field. Rejected fields include all those containing 
a star brighter than 5th magnitude, fields containing cores of 
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Fig. 1. — Distribution as a function of the proper motion for stars found 
with SUPERBLINK. The number density very closely follows a W ~ iT^ law 
(dashed Une, this is not a fit). 



bright globular clusters, parts of M31, and of a few other ex- 
tended and saturated objects. 

Stars identified by SUPERBLINK with proper motions in 
the 0.15" yr"' < n < 2.0" yr"' range and confirmed by visual 
inspection, are found to disproportionately consist of slower- 
moving objects. As the proper motion of a star is inversely 
proportional to the distance, a uniform density distribution of 
stars in the volume around the Sun is expected to result in a 
cumulative distribution inversely proportional to the cube of 
the propermotion. The objects identified with SUPERBLINK 
very closely follow such a A'^ ~ /z"^ relationship, as illustrated 
in Figure 1. The sharp drop in objects above /i = 2.0" yr"' and 
below /i = 0.15" yr"' simply results from the detection limits 
imposed on SUPERBLINK. While the upper limit was firmly 
set into the software, the lower hmit has been set only for the 
purpose of the present catalog. We allowed SUPERBLINK 
to identify stars with proper motions as small as /j, = 0.04" 
yr"' (totaling nearly one million objects), but only those with 
H < 0.15" yr"', considered the most valuable, were examined 
visually and retained for the present analysis. The much more 
numerous slower-moving objects are only now being exam- 
ined, and their publication is planned for a future release. 

In the area analyzed by SUPERBLINK, the software iden- 
tified a total of 56,238 objects with proper motion exceeding 
0.15" yr"'. Among these, a total of 1,159 objects were sub- 
sequently found to be double stars (see §3.2 below), and are 
listed as two distinct objects in the LSPM catalog. This makes 
a total of 57,397 individual high proper motion objects iden- 
tified with SUPERBLINK. The distribution on the sky is dis- 
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Fig. 2. — Top: distribution on tlie north celestial hemisphere of the 57,763 high proper motion stars identified in the Digitized Sky Siirveys with the 
SUPERBLINK software. Bottom: distribution as a function of optical V magnitude. The detection efficiency of SUPERBLINK exceeds 99% in the magnitude 
range 12.0 <V < 19.0. The efficiency drops for brighter (V < 12) stars as the stellar images become saturated on the POSS plates, and at fainter magnitudes 
(V > 19) as one reaches the POSS plate limit. The turnover in the distribution beyond V = 16 is real, and is not a result of a decUning detection rate. It occurs 
because proper motion selected samples survey a limited volume, combined with the fact that field stars also have a turnover in their luminosity function. 



played in Figure 2, along with their distribution as a function 
of optical V magnitude (see §4.4 for a discussion on how V 
magnitudes are derived for SUPERBLINK detections). 

At magnitudes fainter than V = 19, we observe a sharp 
drop in the number of high proper motion stars detected with 
SUPERBLINK. This reflects the limited capabihties of SU- 
PERBLINK to detect stars near the magnitude limit of the 
POSS-I plates. While it is true that the POSS-II plates are 
marginally more sensitive, SUPERBLINK demands a detec- 
tion at both epochs in order to identify the object as moving, 
and thus the detection threshold of SUPERBLINK is deter- 
mined by the sensitivity of the POSS-I plates. 

At the bright end of the distribution {V < 12.0), we observe 
a steady decline in the number of stars detected, falling to 
zero for V < 8. While one naturally expects to find fewer 
high proper motion stars with very bright magnitudes, many 



more high proper motion stars are known with (V < 12.0) 
than have been detected by SUPERBLINK. The lack of stars 
detected at bright magnitudes reflects the inability of SU- 
PERBLINK to deal with stars that have a strong saturated 
core on the photographic POSS-I and POSS-II plates. Tests 
made with fields containing known, bright high proper motion 
stars showed that SUPERBLINK generally did well with stars 
fainter than V = 10.0, correctly identifying them. However, 
tests revealed that SUPERBLINK has much trouble identify- 
ing brighter stars. What generally happens is that the BLINK 
procedure typically fails to determine the centroid of bright 
saturated stars, and is thus unable to calculate a proper mo- 
tion; such objects are simply rejected by the code. In the most 
extreme cases, i.e., for stars brighter than V = 4, a complete 
failure occurs in the SUPER procedure: the code is unable to 
superpose the two 17' x 17' subfields that display extended. 
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saturated patches from the bright star. These fields are ba- 
sically not processed by the code. Unfortunately, this also 
guarantees that any fainter high proper motion star that would 
normally be detected by SUPERBLINK, but happens to be in 
a subfield occupied by a very bright star, will also be missed 
by the code. 

We identified all possible counterparts of the SU- 
PERBLINK objects in the HIPPARCOS and TYCHO-2 cat- 
alogs (see §3.3 below), and adopted for those stars the more 
precise proper motion value from the TYCHO-2 catalog. As a 
result, 91 SUPERBLINK objects that were initially above our 
proper motion threshold were found to have TYCHO-2 proper 
motions below ji = 0.15" yr"'. These rejected objects are not 
considered in the current analysis, and are not counted among 
the 57,397 stars officially identified with SUPERBLINK. 

The 57,397 SUPERBLINK stars form the core of our new 
LSPM catalog, and include thousands of newly identified high 
proper motion stars. 

3.2. Resolved common proper motion doubles 

Common proper motion doubles with separations in the 
range « 1 - 10" are not uncommon in the field. These objects 
usually show up on the POSS plates as elongated, or oddly 
shaped objects, and can be mistaken for distant galaxies or 
short asteroid tracks. Our SUPERBLINK software identifies 
all moving objects, regardless of their shape, and so it picks 
out barely resolved common proper motion doubles just as 
well as single stars. Upon visual inspection, any moving ob- 
ject with an odd shape is flagged for further analysis. 

Objects flagged as possible common proper motion doubles 
are searched for in the 2MASS All Sky Point Source Catalog 
to see whether they are featured as pairs of stars. The res- 
olution of the 2MASS infrared CCD images is significantly 
better than the POSS plates, and pairs of objects with sep- 
arations smaller than 1" are often resolved. We found that 
the vast majority of the stars that were initially flagged by us 
as possible doubles indeed did show up as pairs of resolved 
stars in the 2MASS catalog. All the objects found by SU- 
PERBLINK were eventually searched in the 2MASS catalog, 
in order to determine their infrared (J,H,Ks) magnitudes (see 
§4.2 below). In this process, several more objects identified 
by SUPERBLINK, and not initially flagged by us as candi- 
date doubles, were also found to be resolved into pairs in 
2MASS. The individual components of those close common 
proper motion doubles are included as separate entries in the 
catalog. 

Using all available images from the DSS (POSS-1, POSS-11 
in three colors) plus the 2MASS Quicklook Images obtained 
from the NASA/IPAC Infrared Data Archive ^, we exam- 
ined all the pairs to determine whether they were actual com- 
mon proper motion doubles, or chance alignments. We found 
about equal numbers of each. In areas with significant crowd- 
ing (low-galactic latitude fields) there were an abundance of 
cases in which the high proper motion star happened to be in 
the vicinity of a background source; these were easily filtered 
out by noticing that the background source had not moved be- 
tween the POSS-l/POSS-11 and 2MASS epochs. In most other 
cases, it was clear from the POSS scans and 2MASS images 
that we were dealing with a moving pair There were only a 
few ambiguous cases, for which we conservatively assumed 
the star to be single. 

In the end, we resolved 1,159 SUPERBLINK objects into 

^ http:irsa.ipac.caltech.edu/applications/2MASS/QL/interactive.html 



common proper motion pairs; each component is included in 
the LSPM catalog as a separate entry. However, since it was 
generally not possible to obtain proper motions for each com- 
ponent individually (SUPERBLINK only gave a proper mo- 
tion for the pair) the two stars are listed as having exactly the 
same proper motion. This, of course, is only approximate, 
and one should not conclude that the two stars do not show 
any significant relative proper motion. 

A number of common proper motion doubles were already 
listed as such in the NLTT catalog. In those instances, we have 
tried to assign the correct NLTT numbers for each star of the 
pair To check our assignments, we have first used the coor- 
dinates listed in the NLTT to determine which star of the pair 
was to the north or east of the other. In many cases, the two 
stars were listed in the NLTT as having exactly the same posi- 
tion. In those cases, we looked for notes to the NLTT catalog, 
which usually specified the position angle of the secondary. 

A separate paper (Lepine et al., in preparation) will provide 
a detailed analysis of all the common proper motion doubles 
identified in our survey. 

3.3. Additional Stars from the TYCHO-2 and ASCC-2.5 
catalogs 

Our SUPERBLINK survey of the POSS plates has a bright 
magnitude limit that limited our identification of very bright 
(y < 12) high proper motion stars. In order to build a catalog 
that is the most complete possible, we need to complement the 
SUPERBLINK stars with lists of known, bright high proper 
motion stars. 

The two sources we used to complement our catalog are 
the Tycho-2 Catalogue of the 2.5 Million Brightest Stars 
(TYCHO-2), and the All-sky Compiled Catalogue of 2.5 mil- 
hon s tars (ASCC-2.5). The TYCHO-2 catalogue Wioz etal\ 
l2000h is the product of a re-analysis of data from the ESA 
Hipparcos satellite, and combines space-determined positions 
and proper motions for 2.5 million of the brightest stars in 
the sky (the catalog is complete down to about VV = 12) 
with ground based astrometry from a variety of sources. The 
ASCC-2.5 (Kharchenko 2001) is a catalog largely built from 
the TYCHO and HIPPARCOS catalogs, and providing essen- 
tially similar information on positions and proper motions. 
However, the ASCC-2.5 includes complete data on a number 
of stars whose proper motions and/or photometric data were 
missing in the TYCHO-2 (including stars from the TYCHO- 
2 supplement- 1 , which contains all HIPPARCOS stars not 
Hsted in the TYCHO-2 catalog). The ASCC-2.5 also in- 
cludes astrometric information (including proper motions) on 
an additional number of fainter stars, obtained from various 
ground-based astrometric surveys. The ASCC-2.5 extends the 
TYCHO-2 catalog down to slightly fainter magnitudes. 

The TYCHO-2 catalog was used as our primary source 
of bright, high proper motion stars, while the ASCC-2.5 
was used as a complement to the TYCHO-2. We extracted 
from TYCHO-2 all stars listed with proper motions exceed- 
ing /i = 0.15" yr"' . We found 8,225 objects in the northern sky 
spanning a range in magnitude 2.1 < V < 13.6 (see §4.4 for 
our derivation of V magnitudes from the TYCHO-2 VV and 
Bt magnitudes), with 93% of the stars brighter than V = 12.0 
(Figure 3). Scanning the ASCC-2.5 catalog for additional ob- 
jects, we identified 5,239 stars stars listed with a proper mo- 
tion fj, > 0.15" yr"' that were not Hsted in the TYCHO-2, 
or whose proper motion data were unavailable in TYCHO- 
2. The additional ASCC-2.5 stars spanned a range in mag- 
nitude 0.0 < y < 15.5, but with 95% of the objects fainter 
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Fig. 3. — Distribution as function of magnitude for stars listed in tlie 
TYCHO-2 and ASCC-2.5 catalogs with proper motions larger than 0.15" 
yr"' . The dashed lines show stars listed only in TYCHO-2, while the full line 
shows the full set of TYCHO-2 stars complemented by additional objects 
from the ASCC-2.5 catalog. The complete set is most probably complete 
down to y = 10, but has a sharp drop in completeness fainter than V = 12 at 
which point, fortunately, the detection efficiency of SUPERBLINK reaches 
high levels (see Figure 2). 



than y = 10.0. The vast majority of the V < 10 stars in the 
ASCC-2.5 catalog are listed in the TYCHO-2. 

A comparison of the list of bright TYCHO-2/ASCC-2.5 
high proper motion stars with the list of SUPERBLINK de- 
tections indicated that 4,252 of the TYCHO-2 stars, and 3,603 
of the ASCC-2.5 stars were already in the list of /U > 0.15" 
yr"i stars detected by SUPERBLINK. Among the stars that 
were not found in the SUPERBLINK list, there were signif- 
icant numbers of objects with a TYCHO-2/ASCC-2.5 proper 
motion close to the 0.15" yr~' limit of our initial list of SU- 
PERBLINK detections. We thus surmised that that some 
might have been detected by SUPERBLINK but ranked in a 
lower proper motion range. This is especially true of brighter 
(y < 10) stars that are strongly saturated on the POSS plates 
and thus have larger SUPERBLINK proper motion errors. We 
searched for possible matches in a prehminary list of stars 
found by SUPERBLINK with calculated proper motion 0.10" 
yr"' < /i < 0.15" yr"' . We found matches to an additional 846 
TYCHO-2 and 21 ASCC-2.5 stars. Because our TYCHO-2 
sample contains brighter stars on average than our ASCC-2.5 
sample, it comes as no surprise that most of the additional 
matches were from TYCHO-2 objects (whose proper motion 
errors in SUPERBLINK are larger). 

This left us with 3,127 stars from TYCHO-2 and 1,615 stars 
from ASCC-2.5 that had not been detected by SUPERBLINK. 
Upon visual inspection of DSS/XDSS images centered on 
those objects, we discovered that a significant fraction of them 
do not show any detectable proper motion. This suggests that 
some of the TYCHO-2/ASCC-2.5 stars had their proper mo- 
tions overestimated, which would expliiin why they were not 
identified by SUPERBLINK. 

In order to verify the proper motions quoted in the TYCHO- 
2 and ASCC-2.5 catalogs, we identified counterparts for all 
these objects in the 2MASS All-Sky catalog. We then re- 
calculated their proper motions using the 2MASS positions 



(epochs 1997-2001) and the Hipparcos-based positions from 
the TYCHO-2 and ASCC-2.5 catalogs (epoch 1991.25). The 
2MASS positions for these bright stars are accurate to about 
120 mas, which means that it should be possible to derive 
proper motions to an accuracy of 12 mas yr"' (for a 10 years 
baseline) to 20 mas yr"' (for a 6 years baseline). This as- 
sumes, of course, that the TYCHO-2 and ASCC-2.5 positions 
are significantly more accurate than the 2MASS positions. 

Results showed that a few hundred of the TYCHO-2 stars, 
and over a thousand of the ASCC-2.5 stars had been missed 
by SUPERBLINK for a good reason: their actual proper mo- 
tion is definitely below our adopted threshold of 0.15" yr "'. 
Figure 4 compares the proper motions quoted in the TYCHO- 
2/ASCC-2.5 catalog to the proper motions determined from 
the differences in the 2MASS and TYCHO-2/ASCC-2.5 po- 
sitions. We plot the results separately for 4 groups of objects: 
(1) stars from the TYCHO-2 catalog that were also recov- 
ered by SUPERBLINK (top left), (2) stars not listed in the 
TYCHO-2 catalog but listed in the ASCC-2.5, and that were 
identified by SUPERBLINK (bottom left), (3) stars hsted in 
TYCHO-2 that were not recovered by SUPERBLINK (top 
right), and (4) stars in the ASCC-2.5 but not in TYCHO-2, 
and that were not identified by SUPERBLINK (bottom left). 
First of all, the top left and lower left plots show that there is 
a good correlation between the 2MASS-derived proper mo- 
tions and those quoted in TYCHO-2 and ASCC-2.5, at least 
for stars that had their proper motions confirmed by SU- 
PERBLINK. We find a dispersion of \5mas yr"' in the differ- 
ence between the TYCHO-2 and 2MASS -derived proper mo- 
tions, and ISmas yr"' in the difference between the ASCC-2.5 
and 2MASS-derived proper motions, in good agreement with 
the predicted values (see above). There are very few outliers 
in the distribution, with perhaps a few dozen stars (out of sev- 
eral thousand) whose 2MASS-derived proper motion appears 
to be clearly overestimated, which can be accounted for by an 
occasional, large error in the 2MASS, or TYCHO-2/ASCC- 
2.5, position. 

The upper right and lower right plots, on the other hand, tell 
quite a different story. A significant number of stars are found 
to have 2MASS-derived proper motions around and below 50 
mas yr"' . Given the positional errors on the 2MASS positions, 
these values are consistent with the stars having no detectable 
proper motions. We have visually inspected DSS/XDSS pairs 
of images for over a hundred stars rejected in the procedure, 
and confirmed that indeed none of these stars showed any sig- 
nificant proper motion. The dashed Unes in Figure 4 shows 
where we have set the limits under which a star is consid- 
ered to have no measurable proper motion, in which case the 
quoted TYCHO-2/ASCC-2.5 proper motion is in error. A 
total of 230 TYCHO-2 and 917 ASCC-2.5 presumed high 
proper motion stars were thus determined to be actual low 
proper motion objects. 

Erroneous proper motion entries in the TYCHO-2 catalog 
are certainly cause for concern. We note that most of them 
are associated with stars near the faint end of the TYCHO-2 
catalog (Figure 5); most erroneous entries have V k, 12. We 
plot in Figure 5 the fraction of TYCHO-2 stars with quoted 
proper motion /x > 150mas yr"' that actually made it into the 
LSPM catalog. One can see that fully 20% of V > 12 stars 
were found to be low proper motion objects. 

Furthermore, we could not find any 2MASS counterpart for 
230 of the ASCC-2.5 stars. Visual inspection of DSS/XDSS 
images showed no trace of these stars at the position quoted in 
the ASCC-2.5. Furthermore, no high proper motion star was 
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Fig. 4. — Comparison between the proper motions quoted in tlie TYCHO-2/ASCC-2.5 catalogs (abscissa) and the proper motions calculated from the 
difference between the 2MASS position (epoch 1997-2001) and the TYCHO-2/ASCC-2.5 position ("epoch 1991.25). Top left: 5,098 high proper motion stars 
from TYCHO-2 that have been recovered by SUPERBLINK. Bottom left: 3,624 stars from the ASCC-2.5 (but not listed in TYCHO-2) that have been recovered 
by SUPERBLINK. Top right: 3,127 stars from the TYCHO-2 catalog that have not been recovered by SUPERBLINK (mostly because they are too bright for 
SUPERBLINK to handle). Bottom right: 1,615 stars from the ASCC-2.5 catalog (but not in TYCHO-2) that were not recovered by SUPERBLINK. A significant 
fraction of the TYCHO-2/ASCC-2.5 stars that were not recovered by SUPERBLINK are found to have 2MASS-derived proper motions inconsistent with then- 
quoted values (areas below dashed lines). These have not been included in the LSPM catalog. 



found within 2' of the quoted position. We thus assume these 
entries to be bogus, although we cannot rule out the possibility 
of a very large error (several arcmin) in the quoted ASCC-2.5 
position. 

In addition, there were 66 TYCHO-2 stars and 239 ASCC- 
2.5 stars that are Usted as close visual companions of brighter 
TYCHO-2 objects. These are not resolved on the 2MASS im- 
ages, and thus have no 2MASS catalog counterparts. They 
are, of course, not resolved on the DSS/XDSS images ei- 
ther, and we thus cannot obtain an independent confirma- 
tion of their existence (we also found no mention of them 
in the Luyten catalogs). We refrain from including them in 
the LSPM catalog at this point, while we are still investigat- 
ing their status. We do point out that at least all the primary 
components are in the LSPM, which should make these sec- 
ondaries (and probably many more unsuspected ones) easy to 



recover eventually. 

In the end, we were left with 2,831 stars from TYCHO- 
2, and 229 stars from ASCC-2.5 which are bona fide high 
proper motion stars. Each of these stars has been included in 
the LSPM catalog. The vast majority of those additions are 
stars brighter than V = 10. Their distribution on the celes- 
tial sphere is shown in Figure 6. We find that 271 of the ad- 
ditional TYCHO-2/ASCC-2.5 stars are located in areas that 
were not processed by SUPERBLINK. These areas include 
the north polar cap, and all areas containing very bright stars 
that could not be overlapped properly with the SUPER proce- 
dure. The distribution of those stars as a function of V magni- 
tude is plotted separately in Figure 6. It is obvious that many 
of the rejected areas are coincident with very bright stars, 
as the rejected fields contain a disproportionate number of 
V <5 objects. We finally find 2,775 stars with proper motion 
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Fig. 5. — Top left: distribution as a function of V magnitude of TYCHO-2 stars that have not been included in the LSPM catalog (see Figure 4). Top 
right: the same, for additional stars from the ASCC-2.5 catalog, i.e., stars in the ASCC-2.5 catalog listed with ^ > 0.15" yr"' and which have no TYCHO- 
2 counterparts. The fractional contribution of these misidentified high-/i stars to the full TYCHO-2/ASCC-2.5 samples is shown in the bottom plots. Up to 
rsi 40% of the TYCHO-2 stars listed with jj. > 0.15" do not actually have large proper motions. While a very significant fraction (fs 65%) of the 10 < V < 12 
ASCC-2.5(non-TYCHO-2) stars are bogus, almost all of the V > 12 object are correctly identified (note that most of the latter have also been identified with 
SUPERBLINK). 



// > 0.15" yr' in fields that were analyzedby SUPERBLINK 
but that were missed by the code. Their distribution peaks at 
a magnitude V ~ 8.5, and spans a range in magnitude where 
the efficiency of the SUPERBLINK software is limited be- 
cause of saturation on the POSS plates. 

3.4. Additional stars from the LHS andNLTT catalogs 

A comparison with the LHS and NLTT catalogs reveals a 
small number of stars that are absent from our list of SU- 
PERBLINK detection, and that are too faint to be in the 
TYCHO-2 and ASCC-2.5 catalogs. All the objects were in- 
vestigated individually in order to determine whether they are 
real, and whether they should be added to the LSPM catalog. 

First of all, the LHS includes a list of 13 faint stars with 
proper motions exceeding /i = 2.0" yr"' ; much too fast to have 
been detected with SUPERBLINK. Generally, LHS stars with 



very large proper motions also happen to be relatively bright, 
and are thus also present in the TYCHO-2 catalog. But those 
13 LHS stars are fainter than V = 14, which explains why 
they are not in the TYCHO-2/ASCC-2.5 catalogs either. In 
any case, they were easily re-identified by direct inspection of 
DSS scans, and have been added to the LSPM catalog. 

Secondly, and most importantly, a significant number of 
NLTT stars that are not in the LSPM were actually recovered 
by SUPERBLINK, but found to have proper motions below 
the 0.15" yr~' cutoff of the current version of the LSPM cat- 
alog. There are also 214 bright NLTT stars that are listed in 
TYCHO-2 as /z < 0. 15" yr"' stars. From a search of a prehm- 
inary fist of 0.10" yr"' < fj, < 0.15" yr"' high proper motion 
stars identified with SUPERBLINK, we also recovered an ad- 
ditional 1,204 stars from the NLTT catalog. A small frac- 
tion of these stars were already listed in the NLTT catalog as 
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Fig. 6. — High proper motion stars listed in the Tycho-2/ASCC-2.5 catalogs that were not recovered by SUPERBLINK. Top: distribution on the northern 
celestial hemisphere. Center: distribution as a function of magnitude of the 278 stars that were in areas of the northern sky not processed with SUPERBLINK 
(325 square degrees). Bottom: distribution as a function of magnitude of the 3,787 stars missed by SUPERBLINK because they were too bright, and their images 
saturated on the DSS scans. These additional bright high proper motion stars have all been incorporated into the LSPM catalog. 



having low proper motions, but most were hsted with having 
proper motions 0.18" yr"^ < jj. < 0.25" yr"^ 

Thirdly, we found 161 additional NLTT stars that are ac- 
tually listed in the NLTT as close companions (< 10") of 



brighter NLTT objects. These companions are not resolved on 
the POSS scans, and neither are they resolved on the 2MASS 
images, which explains why they didn't show up in our initial 
search. However, Luyten provided separations and position 
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angles for most common proper motion doubles in the NLTT 
in a comment line, which is available in the electronic version 
of the catalog. Using this information, we rederived the loca- 
tions of those companions (using the revised positions of the 
primaries) and included the companions as separate LSPM 
entries. Exceptionally for those secondaries, the quoted op- 
tical b and r magnitudes (in our LSPM catalog) are directly 
recopied from the NLTT catalog. 

While investigating potential NLTT secondaries that might 
have been missed in our initial search, we also found 39 du- 
plicate entries. These stars are apparently objects whose po- 
sitions and/or proper motions have been remeasured at some 
point, but for which the initial, erroneous entry had been kept 
in the NLTT catalog by mistake. They can be easily identified 
in that while the two stars were supposed to be two objects 
of equal magnitude and colors within 20 - 60" of each other, 
only one object of the pair was found on the DSS scans. The 
fact that the two entries are listed in the NLTT catalog with ex- 
actly the same magnitude and color betrays the fact that they 
are indeed one and the same. 

A more time-consuming job was to investigate the existence 
of the remaining 900 or so NLTT stars that remained unac- 
counted for. Our methodology was simple: we retrieved pairs 
of 17' X 17' DSS/XDSS scans centered on the quoted posi- 
tions of the NLTT stars. The pairs were aligned (shift-rotated) 
using SUPERBLINK subroutines, and examined by eye, by 
blinking them on the computer screen. We searched for the 
presence of a moving object within a 5' x 5' area centered 
on the presumed location of the NLTT star In 395 fields, no 
moving object could be found whatsoever In those instances, 
we must assume that the NLTT entry is bogus. In 7 more 
cases, we did recover a moving object, but its proper motion 
was very clearly smaller than our 0.15" yr~' limit. All of 
these NLTT stars, none of which were included in the LSPM 
catalog, are listed in a separate table in the appendix. 

Finally, we positively identified a total of 492 NLTT stars 
from the DSS/XDSS scans, with proper motions in the 0.15" 
yr"^ < fj, < 0.4" yr"' range (Figure 7). Of those genuine high 
proper motion stars, a total of 174 were found to be located 
in areas that were initially mishandled and rejected by SU- 
PERBLINK, and that were thus not part of the survey. In par- 
ticular, there were 40 stars located north of S = 86.5, the area 
near the north celestial pole that was not properly processed. 

The distribution of stars in areas not analyzed by SU- 
PERBLINK very closely follows the distribution of high 
proper motion stars detected by SUPERBLINK in the rest 
of the sky, if TYCHO-2/ASCC-2.5 stars are excluded. In 
Figure 7, we superpose (dotted line) the distribution of SU- 
PERBLINK objects that do not have a TYCHO-2 or ASCC- 
2.5 counterpart over the distribution of additional NLTT stars 
found in the areas not analyzed by SUPERBLINK. The two 
curves are in very close agreement, with a peak around V = 17. 
(The main difference is that the SUPERBLINK distribution 
extends to slightly fainter magnitudes, which is consistent 
with a high completeness of SUPERBLINK over the NLTT 
beyond V = 19.) This is exactly what is expected for stars 
that are missing because they are outside the SUPERBLINK 
survey areas: they should follow the same general magnitude 
distribution as the stars extracted inside the SUPERBLINK 
survey areas. 

The remaining 318 NLTT stars are in areas that were pro- 
cessed by SUPERBLINK, but that were nevertheless missed 
by the code. Several of these elusive objects were within 1' 
of very bright stars at one of the POSS-I or POSS-II epochs. 



or were close to plate edges, or were coincident with local 
plate defects, making their identification difficult. The distri- 
bution of these stars with V magnitude is also skewed toward 
very faint objects (see Figure 7 bottom panel), with a peak 
at magnitude V = 19.5. This marks the range at which SU- 
PERBLINK is beginning to suffer from incompleteness, as it 
reaches the magnitude limit of the POSS-1 plates. 

3.5. Other additional objects 

Two more stars were included in our catalog that are ob- 
jects with very large proper motions that were discovered in 
the past two years. Both are too faint to have been in the 
TYCHO-2 or ASCC-2.5 catalogs, and because they are very 
recent additions, they are, of course, not in the Luyten cata- 
logs either 

The first star is LSPM 1826H-3014, discovered by 
Lepin e et al.\ J2002h . The star has a proper motion = 2.38" 
yr"' and a magnitude V = 19.4. In our catalog, it bears the 
name LSPM J1826H-3014. The star was actually discovered 
in the course of our own survey, but is regarded as a serendip- 
itous discovery: it was not initially identified as a high proper 
motion star by SUPERBLINK, but rather as a pair of variable 
stars within 2' of each other, which we examined further out 
of curiosity. This makes one seriously consider the possibility 
that there still exist faint stars with very large proper motions 
waiting to be discovered. 

The second star is the extremely h igh proper motion objec t 
SO 025300.5H-165258 discovered bv lTeegarden THH il2003l) . 
The star has a proper motion fi = 5.05" yr"' and a magnitude 
V = 15.4, and is identified as LSPM J0253H-1652 in our cat- 
alog. It is believed to be a very nearby star. Its extremely 
large proper motion is beyond the detection limit (2.0" yr"') 
of the SUPERBLINK. What is interesting is that the survey 
by Teegarden ef fl/. (2003) that led to its discovery was ini- 
tially aimed at the identification of solar system objects, and 
uses a temporal baseline of months to a few years. It thus ap- 
pears that a pair of all sky surveys with a short separation in 
time (e.g. 1-2 years) might well lead the way to locating any 
possible remaining faint stars with proper motions in excess 
of 2.0" yr"i. 

We note that most of the L dwar fs and T dwa rfs discovered 
in recent years (f Kirkpatrick et al. 1999, 200^: "H awlev et all 
I2OOI iCutrief a/.. .2003) very probably have proper motions 
within the range of our catalog; however no effort was made 
to include any high proper motion brown dwarf at this point. 
Proper motions have so far been determined only for a small 
number of L and T dwarfs (Dahn et al. 2002), and additional 
work would be required to obtain accurate proper motions for 
most of them. For now, we have limited ourselves to the very 
few L dwarfs that do show up on the POSS-1 and POSS-II 
plates and that were recovered by SUPERBLINK, although 
we do plan to add high proper motion brown dwarfs to the 
LSPM catalog in the near future. 

Adding up all the stars found by SUPERBLINK, those 
retrieved from the TYCHO-2 and ASCC-2.5 catalog, the 
LHS/NLTT stars missed by our code, and the two additional 
objects discussed in this section, we come to a total and final 
tally of 61,618 stars in the LSPM catalog. 

3.6. Counterparts in the UCAC2 astrometric catalog 

The Second U.S. Naval Observatory CCD Astrograph Cat- 
alog (UCAC2) is the second release in an all-sky astromet- 
ric survey of stars in the magnitude range 7.5 < Rp < 16.0 
JZacharias et a/.ll2003ft . The current version lists 48,330,571 
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Fig. 7. — High proper motion stars listed in the LHS and NLTT catalogs of high proper motion stars, but that were not recovered by SUPERBLINK, and are not 
listed in the Tycho-2 catalog either. Top: distribution on the northern celestial hemisphere (Aitkins projection). Center: distribution as a function of magnitude 
of the 174 stars that were in areas of the northern sky not processed with SUPERBLINK (325 square degrees). The normalized distribution of SUPERBLINK 
stars found in the rest of the sky and that are not in the TYCHO-2/ASCC-2.5 catalog is shown for comparison (dotted line). Bottom: distribution as a function of 
magnitude of the 318 stars missed by SUPERBLINK. These stars have all been added to our catalog. 
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stars in the declination range -90 <Decl.< +50, and gives po- 
sitions with an accuracy of 20-70mas (depending on magni- 
tude). It also provides proper motions for all cataloged stars 
with an accuracy ^ 5mas yr"' . 

We found UCAC2 counterparts for 9,151 of the LSPM 
stars. All the counterparts are south of Decl.=+53.23, re- 
flecting the current, limited sky coverage of the UCAC2 (Fig- 
ure 8). Fewer stars are found below a Decl. of 10 degrees, 
where the UCAC2 apparently has a brighter magnitude limit 
(V < 12). Overall, the UCAC2 lists stars down to a magnitude 
V ~ 16.0 but appears to be significantly incomplete for high 
proper motion stars at all magnitudes and positions. 

While we were initially hoping to use the UCAC2 as a an 
additional source of bright stars with large proper motions for 
the LSPM catalog, we found this to be impractical at this 
point. The main reason is that it appears that the UCAC2 
is plagued with a large number of spurious high proper mo- 
tion entries. We retrieved all the stars in the UCAC2 that are 
in the northern sky and listed as having a proper motion ex- 
ceeding 0.15" yr"^ Apart from the 9,168 objects also listed 
in our LSPM catalog, we found an additional 7,370 entries 
(see Figure 6) with a very non-uniform distribution. The vast 
majority of the additional entries are located at low galac- 
tic latitudes, and they have cataloged proper motions in the 
range 0.15" yr"' < < 0.25" yr"'. Our examination of sev- 
eral DSS/XDSS scans centered on the presumed location of 
those stars failed to reveal any of them as high proper motion 
objects. 

Good UCAC2 counterparts do however prove useful as a 
means to estimate the astrometric accuracy of the LSPM from 
an independent source. A comparison of the UCAC2 and 
LSPM positions and proper motions for the stars common to 
both catalogs is presented in §4.3 and §4.4 below. 

4. PHOTOMETRY 

4.1. Optical magnitudes from TYCHO-2 and ASCC-2.5: B 

and V 

Optical Bt and V7 magnitudes are obtained from TYCHO- 
2 catalog counterparts (see §3.3 above). Photometric errors 
are 0.013 mag for Vt < 9, and 0.1 mag for 9 < Vt < 12. 
A simple conversion transforms these into Johnson B and V 
magnitudes: 

B = Bt-0.24(Bt-Vt), (1) 
V = Vt-0.09(Bt-Vt), (2) 

following the prescription in the introduction to the HIPPAR- 
COS and TYCHO catalogs. 

The ASCC-2.5 catalog provides both B and V magnitudes 
(converted from Vt and Bj), and can also be used as a rela- 
tively reliable source of optical magnitudes for bright stars. In 
particular, we are using it to obtain magnitudes of bright stars 
that are not listed in the TYCHO-2 catalog. For stars fainter 
than about V = 12.0, ASCC-2.5 magnitudes are derived from a 
variety of sources, and may not be as accurate as the TYCHO- 
2 magnitudes, but since they were obtained from photoelectric 
or CCD measurements, they should be relatively reliable. 

We have gathered B and V magnitudes from TYCHO- 
2/ASCC-2.5 for a total of 11,719 LSPM stars. The frac- 
tion of stars with TYCHO-2/ASCC-2.5 optical photometry 
is plotted in Figure 9 as a function of V. It shows that the 
LSPM contains reliable optical photometry for essentially all 
stars brighter than V = 12.0. This is fortunate, because these 
are stars for which photographic magnitudes are subject to 



large errors, because of saturation on the POSS plates. For 
fainter stars, and especially those with V > 14.0, we do need 
to rely mainly on photographic magnitudes to cover the opti- 
cal regime: the only existing all-sky catalogs of faint optical 
stars are based on photographic plate material. 

4.2. USNO-Bl.O photographic magnitudes: Bj, Rp, 1^ 

The USNO-Bl.O Catalog (i Monet ef a/.l I200I is an afl- 
sky catalog made from scans of several photographic sky 
survey, including the POSS-1 and POSS-II. Astronomical 
objects have been identified using the PMM scanning ma- 
chine. The catalog gives positions, proper motions, photo- 
graphic magnitudes in 5 passbands, and star/galaxy estimators 
for 1,042,618,261 objects. While the USNO-Bl.O provides 
proper motions for all objects detected in the POSS plates, it is 
not a reliable source for stars with large proper motions. The 
main difficulty with the USN O-B 1 .0 is the exceedingly large 
number of spurious entries ( Gouldll2003bh : at high galactic 
latitudes, up to 99% of objects listed with fx > 0.18" yr"' are 
not real. The catalog also suffers from serious incomplete- 
ness for high proper motion stars at low galactic latitude (up 
to 30%), as estimated from its recovery of NLTT stars. 

The USNO-B 1 .0 is however an extremely valuable comple- 
ment to the LSPM catalog because it provides reasonably ac- 
curate photographic Bj, Rp, and In magnitudes (respectively 
llla-J, llla-F, and IV-N) derived from the POSS-11 scans. Be- 
cause the USNO-Bl.O is based on some of the same plate ma- 
terial as the DSS, it also provides a very useful check for our 
SUPERBLINK detections. 

We have succeeded in finding USNO-Bl.O counterparts for 
60,396 of our LSPM stars. Searching the USNO-Bl.O for 
high proper motion objects however turned out to be a dif- 
ficult and time-consuming problem. Cross-correlation of the 
two catalogs yielded only « 75% of unambiguous matches. 
The large number of ambiguous cases fell into three broad 
classes: (1) stars with erroneous USNO-Bl.O proper motions, 
(2) moving stars not identified as such in the USNO-Bl.O and 
listed as separate stars, one for each epoch in which the stars 
was detected, and (3) confusion with background stars at the 
detection epoch. 

One common problem was USNO-Bl.O entries with large 
errors in their quoted proper motions {ij.usno-bi.o)- Since their 
quoted RAJ2000 and DEJ2000 are calculated by extrapolat- 
ing the position from the mean epoch of observation with 
their estimated proper motions, some stars have quoted posi- 
tions incorrect by up to several seconds of arc. In most cases, 
we were able to recover the star by extrapolating back to the 
mean epoch of observations and recalculating the RAJ2000 
and DEJ2000 positions using the proper motion determined 
by SUPERBLINK. 

Another major source of complication was high proper mo- 
tion stars listed as two or more distinct entries in the USNO- 
Bl.O. Each entry typically corresponds to a detection of the 
star in a distinct photographic survey. This problem partic- 
ularly affected stars with larger proper motions. A typical 
example is a moving star listed as three separate entries, one 
with the position of the star at the epoch of the POSS-1 survey, 
one with the position of the star at the epoch of the POSS-II 
red and blue surveys, and one with the position of the star at 
the epoch of the POSS-11 near infrared survey. The confu- 
sion was such that most of these cases had to be dealt with 
individually. 

Additional complications occurred because of confusion 
with background sources. This problem was common espe- 
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Fig. 8. — High proper stars in the UCAC2 astrometric catalog. Top: positions of 9,151 stars with jj, > 0.15" yr ' that are also listed in our LSPM catalog. 
Bottom: positions of 7,370 stars listed in the UCAC2 as having proper motions fj, > 0.15" yr~^ but that are not in the LSPM; these appear to be spurious entries. 



cially in low galactic latitude fields, where crowding is signif- 
icant. All these cases had to be examined one by one. Overall, 
we had to visually inspect w 12,000 LSPM objects in order 
to determine their correct USNO-B 1 .0 counterpart. Again we 
made use of our interactive software, this time overlaying the 
USNO-B 1.0 catalog over our SUPERBLINK finder charts. 
Whenever an LSPM star appeared as two or more separate 
USNO-B 1.0 entries, a choice had to be made as to which one 
should be used as the "official" counterpart, in order to keep 
the LSPM catalog simple. In general, we picked the entry 
having the most coinplete photometric data, or in some cases 
the one least hkely to be contaminated by blending with back- 
ground sources. 

Not all USNO-B 1.0 entries have magnitude information in 
all three bands. Whenever possible, we tried to combine mag- 
nitude data if a USNO-B 1 .0 star appeared as more than one 
entry. For example, if one high proper motion star was listed 
as two distinct USNO-B 1 entries, one giving only Bj and Rp 
magnitudes, the other giving only an In magnitude, we would 



combine the information to obtain complete BjRfIn photom- 
etry. For practical purposes, however, we hst the counterpart 
ID only for one of the two entries. As a result, LSPM mag- 
nitudes are often more complete than magnitudes extracted 
from the USNO-B 1 catalog for the listed counterpart. 

Finally, we note that no USNO-B 1 counterparts could be 
found for a total of 1,580 LSPM stars. The majority of these 
are in close proper motion double systems or in very crowded 
fields, and are simply not resolved in the USNO-B 1.0, al- 
though most of them are resolved in 2MASS. We note that 
267 of the LSPM stars have neither a 2MASS nor a USNO-B 1 
identifier; a majority of these are faint proper motion compan- 
ions that are not resolved in the USNO-B 1 and are too faint 
to have been detected by 2MASS. However, 229 of them are 
NLTT stars, and the others all clearly show up on DSS/XDSS 
scans and/or 2MASS images. 

4.3. 2MASS infrared magnitudes: J, H, and 
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Fig. 9. — Fraction of LSPM catalog stars with B and V magnitudes ob- 
tained from the TYCHO-2 and ASCC-2.5, plotted as a function of the visual 
magnitude. This indicates that we have reliable optical photometry for essen- 
tially all stars brighter than V = 12.0. For fainter objects, we have to rely on 
photographic magnitudes to estimate V. 



The 2MASS All-Sky Point Source Catalog is an all-sky 
catalog of sources detected in the Two Micron All-Sky Sur- 
vey (|Cutri et al. 2003). The catalog covers the whole sky and 
is complete down to 7 ~ 16.5. Infrared J, H, and mag- 
nitudes are provided, and are accurate to 0.02 mag down to 
15th magnitude. Positions are given for the epoch of observa- 
tion (1997-2001) and are accurate to 70-80 mas for the fainter 
sources {J > 9), and 120 mas for the brighter ones. 

The vast majority of the LSPM catalog stars are found to 
have counterparts in the 2MASS All-Sky Point Source Cata- 
log: we have reliable matches for 59,684 of our high proper 
motion stars. Finding 2MASS counterparts was straight- 
forward for « 90% of the LSPM objects, and was ob- 
tained by a simple cross-correlation of the SUPERB LINK and 
TYCHO-2/ASCC-2.5 high proper motion star positions with 
the 2MASS source positions. To make the search more ef- 
fective, all positions were locally extrapolated to the epoch 
of the 2MASS observations. Most 2MASS counterparts were 
found within 1" of their predicted position. There were multi- 
ple possible matches for « 3,000 of the LSPM objects, most 
of them in very crowded, low galactic latitude fields. We used 
again a software package developed by SL, which overlays 
2MASS catalog entries on the 3.25' x 3.25' charts generated 
by SUPERBLINK. All matches were then made interactively, 
by direct visual inspection. All common proper motion dou- 
bles (see §3.2 above) were also examined and their 2MASS 
identification verified with the same software. No 2MASS 
counterparts were found for 2,292 LSPM stars; in the major- 
ity of cases, these are simply too faint in the infrared to have 
been detected in the 2MASS survey. 

Because of its relatively high astrometric accuracy, we have 
adopted the 2MASS catalog as the primary source of posi- 
tional information for stars that are not listed in the TYCHO- 
2 catalog (see §5.3 below). We believe the 2MASS is the 
best possible choice to pinpoint the positions of those stars 
because: 1) it contains the majority of the LSPM objects, 2) 
its positions are given in the ICRS system and are reason- 
ably accurate, and 3) the 2MASS observation epochs are very 



close to 2000.0. The advantage of having all the 2MASS po- 
sitions to within a few years of the 2000.0 epoch is that it min- 
imizes positional errors introduced by the errors in the proper 
motions, when one extrapolates the position of a high proper 
motion star to the 2000.0 epoch. 

4.4. Estimated V magnitudes and V-J colors 

We have attempted to place all our stars on a simple, uni- 
form magnitude/color system, that includes both an optical 
and an infrared magnitude. We do have uniform, reliable mea- 
surements of infrared magnitudes (from 2MASS) for the vast 
majority of LSPM stars. However, some of our objects do 
not have 2MASS counterparts. Furthermore, we are in a sit- 
uation where our optical magnitudes are in two different sys- 
tems. On the one hand, we have V and B magnitudes accurate 
to 0.1 mag, but only for a small fraction of our stars (those 
with TYCHO-2 counterparts). On the other hand, we have 
much less reliable photographic magnitudes (Bj, Rf, 1^), ac- 
curate to wO. 3-0.5 mags (Monet ef a/. 2003), but most likely 
affected by systematics errors (Sesar et al. 2004). Additional 
complications include the fact that one or more of the photo- 
graphic magnitudes ar e sometimes missing. It is claimed by 
ISalim & Gould' (2003') that photographic R magnitudes from 
the USNO-A2.0 are accurate to 0.25 mag, which is signifi- 
cantly better than USNO-B 1 .0. On the other hand, the USNO- 
Bl.O catalog is more complete than the USNO-A2.0, partic- 
ularly for faint {V > 19) stars, and also at low Galactic lati- 
tudes. Nevertheless, we are currently trying to find USNO- 
A2.0 counterparts of LSPM stars, and USNO-A2.0 magni- 
tudes will be included in future versions of the catalog. 

In any case, it is desirable to provide an immediate means to 
classify the stars in our catalog according to color and magni- 
tude (even roughly). The general idea is thus to get estimates 
of the optical V magnitude and of the optical-to-infrared V -7 
color /or all the stars in the catalog. We already have reliable 
V magnitudes for LSPM stars with TYCHO-2 (see equation 
2) or ASCC-2.5 counterparts, and J magnitudes for all stars 
with a 2MASS counterpart. What we need is a transforma- 
tion system to obtain estimates of V and V-J using the pho- 
tographic magnitudes from the USNO-B 1.0 catalog. While 
one might be tempted to use the excellent 2MASS infrared 
magnitudes to obtain estimates of optical V, this cannot be 
done reliably. The reason is that M dwarfs, which constitute 
the vast majority of the LSPM stars with no TYCHO-2 coun- 
terparts, are largely degenerate in their infrared colors: all M 
dwarfs (except the very coolest) have J-Ks ^ 0.7 ± 0.2 for 

3 <y-y<6. 

The V band is located halfway between the photographic 
Bj and Rf bands. We estimate V using: 

V = Bj-Q.A6{Bj-Rf), (3) 

a relationship which is verified for all TYCHO-2 stars with 
USNO-B 1.0 counterparts. For stars with B j, but no Rp mag- 
nitudes, we find it is possible to estimate V from the following 
set of transformations: 

y = By-0.23(By-7)-0.10 [By-7<4] 

V = Bj-Q.Q5{Bj-J)-Q.12 [Bj-J>A] (4) 

Likewise, for stars with Rf, but no Bj, we can estimate V 
using: 

V =Rf + 0.6(Rf-J)-0.\0 [Rf-J<2] 

V = Rf+\.\Q [Rf-J>2] (5) 
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The V magnitudes estimated from the relationships given 
above are generally accurate to about ±0.5 mag. Follow- 
ing these simple transformations, we calculate V magnitude 
estimates for 61,550 LSPM stars. Of the 427 LSPM stars 
for which we do not provide a V magnitude estimates, 355 
are close common proper motion doubles that are resolved in 
2MASS but not on the POSS plates (see §3.2). At this point, 
we refrain from trying to obtain a V magnitudes using only 
J, H, and Kg. The remaining 72 stars with no V estimates 
are stars that are not in the 2MASS of USNO-Bl.O catalogs, 
and for which we only have Rf magnitudes estimates from 
SUPERBLINK. 

Prospective catalog users should be warned that these V 
magnitude estimates are generally not very accurate, and may 
be subject to systematic errors and other effects. The LSPM 
catalog V magnitudes should only be trusted for stars brighter 
than V = 12.0, whose V are from the TYCHO-2 catalog. At 
fainter magnitudes, there may be errors of 0.5 mag or larger. 

Since the majority of LSPM objects do have 2MASS coun- 
terparts, calculations of V-J colors are straightforward. For 
stars that do not have 2MASS counterparts, we use the photo- 
graphic Iff magnitudes and use the following transformation: 

y-y= i.3(y-/A,)+o.3. (6) 

From the 2MASS counterparts and the transformation above, 
we obtain V-J colors for all but 814 entries in the LSPM 
catalog. These include the 427 stars for which we have no V 
magnitude estimates (see above), and 387 stars that are not in 
the 2MASS catalog and for which we do have 1^ magnitudes. 
Note that the V-J colors are only as accurate as the V mag- 
nitudes are. Since stars fainter than V = 12.0 have errors of up 
to 0.5 mag or even larger, the V-J color estimates should be 
used with extreme caution. 

We emphasize again that the primary goal of the LSPM 
catalog is to provide the most complete list possible of high 
proper motion stars, and is not intended to be a photometric 
catalog. The photometry that we do provide for LSPM stars 
should be regarded as very preliminary, and is given only as a 
help in identifying interesting classes of objects for follow-up 
observations. Future efforts will be devoted to obtaining more 
accurate optical magnitude estimates for all LSPM stars. 

5. ASTROMETRY 

5.1. Conversion to absolute proper motions 

The SUPERBLINK software was largely designed to 
achieve the highest possible recovery rate for high proper mo- 
tion stars on photographic plates. A such, it was optimized 
for raw detection, and not for accurate astrometric measure- 
ments of detected objects. The main caveat is that proper 
motions are calculated relative to local background sources. 
This means, typically, all objects within « 4' of the moving 
target. Because of this, the proper motions calculated by SU- 
PERBLINK are local, relative proper motions. These are usu- 
ally offset by up to several milliarcseconds (mas) per year rel- 
ative to absolute proper motions, which are proper motions 
measured in a fixed reference frame (defined e.g. by the posi- 
tions of distant quasars, such as for the ICRS reference frame). 
It must be realized that most "background" sources used by 
SUPERBLINK as a local reference system are Galactic ob- 
jects, and they all have significant proper motions at the mas 
level. The local frames used by SUPERBLINK are thus mov- 
ing frames, and this potentially introduces both random and 
systematic offsets in the SUPERBLINK proper motions. The 
random offsets arise because of the hmited number of stars 



that locally define the frame, if these stars are all moving in 
random directions, then their mean proper motion will gen- 
erally not add-up to zero; it will however converge to zero if 
the number of reference stars is large enough. These random 
errors affect most the fields at high Galactic latitudes, where 
the object density is low and the local SUPERBLINK refer- 
ence frames are defined by very few stars (sometime < 100). 
In any case, field background stars generally have random 
proper motions smaller than 10 mas yr"\ which means that 
local random offsets will be less than 1 mas yr~^ in frames 
defined by at least 100 stars. Systematic offsets, however, are 
more of a problem. If all local slars participate in some local 
bulk motion, then the local frame used by SUPERBLINK will 
definitely be moving at the bulk motion rate, no matter how 
many stars define the frame. 

Fortunately, we do have a means to estimate some of the 
systemic motions of the background stars, and correct for 
them in order to obtain absolute proper motions. To that 
purpose, we can use all TYCHO-2 stars that have also been 
measured with SUPERBLINK, and compare their absolute 
and relative proper motions. The random errors on the SU- 
PERBLINK proper motion are on the order of, or larger than 
the local systemic motions, but we can average out the resid- 
uals over appropriate- sized areas, and calculate zonal correc- 
tions. 

It is true that the averaging procedure may even out some of 
the local fluctuations. However, background Galactic stars are 
expected to display mainly global patterns of systemic mo- 
tions. The main sources are the rotation of the Galaxy, the 
systemic motion of the local standard of rest (LSR) relative 
to other Galactic stellar populations (old disk, halo), and the 
motion of the Sun within the LSR. The resulting systemic ab- 
solute proper motions of the background stars are dependent 
on their position on the sky, but on a global scale. For in- 
stance, the systemic drift of old disk and halo stars is largest 
in a direction perpendicular to the Galactic rotation (toward 
the Galactic pole), and slowly decreases as one looks more 
toward the direction of rotation. 

The current version of the LSPM catalog only lists a few 
thousand stars (with n > 0.15" yr~') for which we have both 
TYCHO-2 and SUPERBLINK proper motions, leaving very 
few objects to calculate zonal corrections on a scale of less 
than a few tens of degrees. This tends to make the map too 
coarse, and local values too inaccurate. But the complete SU- 
PERBLINK database actually includes detections from the 
DSS of stars with proper motions down to 0.04" yr~K While 
most of those lower proper motion detections are still being 
processed and analyzed, we have recently compiled a prelim- 
inary list of objects, which include over 30,000 TYCHO-2 
stars with proper motions 0.05" yr~' < p< 0.15" yr~\ all of 
which have a relative proper motion value measured by SU- 
PERBLINK. 

Zonal corrections have thus been calculated using a list of 
33,312 stars with magnitudes 10 < V < 13 and proper mo- 
tions 0.05" yr"' < p < 0.50" yr"'. For each position of the 
celestial sphere, the zonal correction is calculated from the 
mean of the offsets between the relative (SUPERBLINK) and 
absolute (TYCHO-2) proper motions of all stars within a ra- 
dius of 7 degrees. Every position on the sky uses « 290 stars 
on average. All offsets are calculated in the local plane of the 
specified location. Outliers, with offsets more than 3-sigma 
away from the mean, are removed from the final calculation. 

Depending on the position on the sky, the mean offsets 
vary from -9.1 to -1-12.2 mas yr"^ in pJiA, and from -0.3 to 
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Fig. 10. — Local offsets between SUPERBLINK relative proper motions and TYCHO-2 absolute proper motions. The offsets are calculated using 33,300 
TYCHO-2 stars with proper motions 0.04" yr"' < fi < 0.50" yr"' whose relative proper motions have been independently calculated with SUPERBLINK. The 
three bottom panels shows local differences between the relative and absolute proper motions for stars located within 7 degrees of the specified location. Proper 
motion differences are calculated in the local plane of the sky. The crosshairs mark the mean value of the offset, while the ellipse shows the mean standard 
deviation (Icr). The top panel plots proper motion difference vectors as a function of position on the sky, in galactic coordinates. The apex (0) and antapex 
(0) of the Sun's motion are noted. These offsets effectively map out the local mean proper motion of background field stars in the TYCHO-2 (ICRS) reference 
system. Offsets are largest at high galactic latitude, where there is a significant drifting motion of old disk and halo star's relative to the local standard of rest. 



+13.9 mas yr~ in ^DE. A map of the zonal corrections for 
the northern sky is shown in Figure 10, where it is plotted 
in Galactic coordinates. The local distribution of offsets be- 
tween the relative and absolute proper motions is also shown 
for three positions on the sky. The number of stars used in 
calculating the local offset is noted. 

From these zonal corrections, absolute proper motions are 
calculated from the relative SUPERBLINK proper motions. 
We thus obtain absolute proper motions for all LSPM stars. 
In the LSPM catalog, we list both the relative proper motion 
determined with SUPERBLINK, and the absolute proper mo- 
tion obtained after applying the zonal corrections. 

In effect. Figure 10 plots the local mean value of the abso- 
lute proper motion of background stars (N.B.: the zonal cor- 
rection vectors plotted in Figure 10 all point in the direction 
opposite to the mean relative proper motion). Ou r plot should 
be compared to Figure 13 in' Munn et al.\ J200 4V which plots 
local mean values of the proper motion of background stars, 
calculated by combining astrometry from the USNO-B 1 .0 and 
SDSS Data Release 1. Exactly the same pattern is displayed 
in both figures. Our Figure 10, however, covers a much large 



area on the sky (20,000 square degrees, compared to the 3,000 
square degrees of the Munn et al. survey), and conveys a more 
global picture of the kinematics of the local Galactic stars. 

The patterns observed in Figure 10 are most probably a 
combination of three effects. The dominant pattern appears to 
be the drift of the LSR relative to the stars in the old disk and 
halo. The mean motion of the standard of rest relative to the 
centroid of the velocity distribution of all local Galactic stars 
is pointing in the direction of galactic rotation (l=90,b=0), and 
there should thus be a general drift of the background stars in 
a direction opposite to the Galactic rotation (l=270,b=0). This 
is largely what is observed in Figure 10, although the relative 
proper motions diverge from a point that does not exactly co- 
incide with (l=90,b=0). But other systematic motions should 
be producing drifting motions in the background stars. One 
is the motion of the Sun itself relative to the stars in the lo- 
cal standard of rest. It should be producing a general drift- 
ing motion pointing away from the apex of the Sun's motion. 
For comparison, we plot in Figure 10 the positions of the So- 
lar apex/antapex, as determined by F ehrenbach et al, ( 2001i) . 
From Figure 10, it is unclear how much weight this effect has. 
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Another effect is the rotation of the Galaxy as a whole, which 
should result in a rotating motion around the north and south 
Galactic poles. 

A detailed model that would account for all those possible 
effect would be required for a more complete interpretation 
of Figure 10. We note that the information summarized in 
our Figure 10 can all be recovered (possibly with even greater 
detail) from a simple subtraction of the relative and absolute 
proper motions, which are listed in separate columns in the 
LSPM catalog. 

5.2. Accuracy of LSPM proper motions 

A check on the accuracy of our proper motions is obtained 
from LSPM stars that have counterparts in the UCAC2 cata- 
log. We separate these objects into two groups: stars that have 
their proper motions directly transcribed from the TYCHO-2, 
and stars that have their proper motions determined by SU- 
PERBLINK. 

The first group includes 4,181 TYCHO-2 stars (all listed in 
the LSPM catalog) that also have UCAC2 counterparts. As 
expected, proper motions from the TYCHO-2 catalog are in 
very close agreement with the UCAC2 proper motions. The 
dispersion in the difference is 2.8 mas yr"' in R.A. and 2.6 
mas yr"' in Decl. This is comparable to the quoted rms errors 
on the UCAC2 proper motions. 

The second group comprises 4,380 stars from the LSPM 
in the magnitude range 12 < V < 16 that have UCAC2, but 
no TYCHO-2, counterparts. Overall, the difference between 
SUPERBLINK and UCAC2 proper motions has a dispersion 
of [7.5,6.7] mas yr"' in RA and Decl. after removal of 3cr 
outliers (Figure 11). There is also a small offset of [0.9,-0.9] 
mas yr"', which might indicate a problem in the zonal cor- 
rections procedure. Our zonal corrections are based on mo- 
tions from relatively bright (TYCHO-2) stars, while the local 
frames used by SUPERBLINK are largely defined with re- 
spect to the more numerous fainter stars. Perhaps the zonal 
corrections are dependent both position of the sky and mag- 
nitude. If bright stars are on average closer to us than fainter 
background stars, then there could very well be a difference 
in their systemic proper motions. In any case, since the nom- 
inal errors on the UCAC2 proper motions are very small (1-3 
mas yr"'), the measured «7 mas yr"' dispersion is a good es- 
timate of the SUPERBLINK astrometric errors on the proper 
motions. 

The 3(7 outliers comprise 6% of the stars in the sample, and 
2.5% of the objects are beyond the 6(j limit. This means there 
is an extended tail to the distribution. Indeed, 90 stars have 
a difference in proper motion > 100 mas yr"'. Whether the 
large difference arises from a faulty LSPM or UCAC2 proper 
motion remains to be determined, although we do suspect that 
in a significant number of cases it is the UCAC2 proper mo- 
tion that is in error 

The SUPERBLINK proper motion error appears to be inde- 
pendent of magnitude for stars fainter than V = 1 1 (Figure 12). 
The proper motions of fainter stars are perhaps marginally 
better, and we measure a dispersion of [7.3,6.3] mas yr"' at 
y > 15. Astrometric errors increase significantly for brighter 
(y < 11) stars, as expected from the fact that these are sat- 
urated on the DSS scans. This, however, is of little conse- 
quence for the proper motions quoted in the LSPM catalog 
since we use the more accurate TYCHO-2 astrometry for the 
vast majority of the LSPM stars with y < 1 1 . The accuracy 
of the SUPERBLINK astrometry is also largely independent 
of the proper motion (Figure 13). 



The proper motion errors from SUPERBLINK are rela- 
tively large r than those quo t ed for the revised NLTT catalog 
(rNLTT) of ISalim & Gouldl j2003h . which are claimed to be 
~ 5.5mas yr"' in both RA and Decl. It is possible that SU- 
PERBLINK errors are slightly larger because SUPERBLINK 
uses photographic plate material for both its first and second 
epoch, while Salim & Gould (2003) used data from photo- 
graphic plates only for their first epoch (USNO-A catalog, 
based on POSS-1) while they used the 2MASS Second Incre- 
mental Release as their second epoch. It is also possible that 
the larger errors in the SUPERBLINK proper motions arise 
from small-scale fluctuations in the systemic motions of the 
background stars, which could only be corrected by a higher- 
resolution map of zonal corrections. A more likely possibility 
is that the SUPERBLINK proper motions are affected by sys- 
tematic errors introduced by astrometric magnitude equations 
(see §5.6 below). 

Readers interested in having more accurate proper mo- 
tions may want to ch eck if their star is in the rNLTT catalog 
(Salim & Gould 2003). The recovery of rNLTT proper mo- 
tions is straightforward, since both the rNLTT and LSPM cat- 
alogs provide NLTT identification numbers. One limitation is 
that the rNLTT only has data for 15,899 of the northern NLTT 
stars, or roughly a quarter of the LSPM stars. 

Clearly, there is still room for improvement, and future ef- 
forts will be devoted to obtain more accurate proper motion 
measurements, which will be included in future versions of 
the LSPM catalog. The possibility of obtaining much more 
accurate proper motion measurements using data from CCD- 
based surveys, such as the Sloan Digital Sky Surve y, was 
demonstrated recently (Gould & Kollmeier 2004; Munn etalX 
I2OOI . Careful astrometric calibration using local quasars can 
yield proper motions with an accuracy < 4mas yr"' . 

The goal of the current version of the LSPM is to provide 
the most complete list of objects possible, with reasonable 
astrometric accuracy. Though significant improvements of the 
proper motion errors are possible, at least for a fraction of the 
LSPM stars, this will require substantial efforts, which are 
beyond the scope of this paper In any case, our positions and 
proper motions are accurate enough to provide a solid starting 
point for future improvements. 

5.3. Accuracy of LSPM 2000.0 positions 

The brighter LSPM stars have their 2000.0 positions ex- 
trapolated from the 1991.25 positions and proper motions 
of their TYCHO-2 counterparts. Fainter stars with no 
TYCHO-2 counterparts have their 2000.0 positions extrap- 
olated from the positions of their 2MASS counterparts and 
their SUPERBLlNK-derived absolute proper motions. 

We estimate the positional accuracy of the fainter (V > 12) 
LSPM objects by comparing the SUPERBLlNK-derived po- 
sitions to those of the UCAC2 catalog (see §3.6 above), 
for those stars that have UCAC2 counterparts (Figure 
14). The difference in position has a calculated dispersion 
[91,88](mas) in [RA.DE]. Note that UCAC2 has a reported 
astrometric precision of 20-70 mas in that range of magni- 
tudes. There is also an offset [-7.6,6.8](mas), which is small 
compared to the magnitude of the dispersion, but is statis- 
tically significant. In principle, the offset could be due to 
small systematic errors in the SUPERBLINK proper motions, 
which are used to extrapolate the 2MASS positions to the 
2000.0 epoch. However, the [0.9,-0.9] mas yr"' system- 
atic offset between the SUPERBLINK and UCAC2 proper 
motions is not consistent with the [-7.6,6.8](mas) positional 
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Fig. 1 1 . — Difference between the SUPERBLINK and UCAC2 absolute proper motions. The zonal-correction has been applied to the SUPERBLINK proper 
motions. Only stars that have no TYCHO-2 counterparts are shown. There is a dispersion of 7.9 mas yr"' in fiRA and 7.1 mas yr^' in fiDE. This provides an 
estimate of the rms errors in the SUPERBLINK proper motions. 



offset. In any case, we find that our positional errors are 
marginally larger than those reported for the entire 2MASS 
catalog, which have a dispersion w 80 mas relative to the 
UCAC2 catalog. 

Most of the brighter LSPM stars (V < 12) have their 
2000.0 positions extrapolated from the TYCHO-2/ASCC-2.5 
positions (epoch 1991.25), using the TYCHO-2/ASCC-2.5 
proper motions. A comparison of the 2000.0 positions with 
those given in the UCAC2 catalog shows a dispersion of 
[81,71](mas) in [RA,DE]. One must note that the positional 
errors in the 2MASS catalog are significantly larger for bright 
stars (> 120 mas forK, < 8). The TYCHO-2 positions there- 
fore remain the best choice at this time. 

5.4. Comparison with the NLTT catalog 

The LSPM catalog provides new estimates of the positions 
and proper motions of all northern stars in the NLTT cata- 
log. As demonstrated by SaHm & Gould (2003) . the NLTT 
positions are accurate to no better than a few arcseconds, with 



some stars having positional errors of up to a few arcminutes. 
The LSPM positions, which are accurate to within one arc- 
second (see above) are a significant improvement. Improved 
position s for 3 1,000 NLTT stars are also available from the 
rNLTT dSahm & Gould 2003) . but only for a little more than 
half the NLTT stars in the northern sky, while our LSPM is 
complete for northern NLTT stars. 

The difference between the NLTT and LSPM proper mo- 
tions shows a dispersion [20.3, 18.7] mas yr"' in [pRA^jiDE} 
(after removal of 3(T outliers), with an offset [1.1,4.4] mas 
yr"' (Figure 15). This estimate of the NLTT errors is 
comparable to the value of cr ~ 20 mas yr"' estimated in 
Salim & Gould (2003), which was based on a comparison be- 
tween NLTT and rNLTT proper motions. The [2.3,6.6] mas 
yr"' offset is also the result of NLTT proper motions being 
relative. Indeed, if we compare NLTT proper motions to the 
relative proper motions measured with SUPERBLINK, the 
offset is reduced to [0.2,-1.6] mas yr"'. 

The accuracy of NLTT proper motions is naturally expected 
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Fig. 12. — Dispersion in the difference between SUPERBLINK and 
UCAC2 proper motions, as a function of magnitude. The continuous line 
shows the mean of the difference, and the dotted lines show the 1 cr disper- 
sion. All stars with SUPERBLINK and UCAC2 proper motions have been 
included here, including those which also have TYCHO-2 counterparts (to 
increase the sample of objects at bright magnitudes). For stars brighter than 
V = ll.O the accuracy of the SUPERBLINK proper motions degrades as the 
stars gets brighter, but the mean errors are uniform in the 11.0 < V < 17.0 
range. 



to be larger, since the second epoch of Luyten's survey was 
only « 15 years after the POSS-1 epoch, while the second 
epoch of the DSS (POSS-11) is 40 years after that of the POSS- 
I. The fact that the NLTT proper motion errors are approx- 
imately three times as large as the SUPERBLINK errors is 
consistent with the difference in the temporal baseline. In any 
case, the LSPM proper motions are a significant improvement 
over NLTT proper motions. Not only is the accuracy better by 
a factor three, but the LSPM proper motions are absolute, in- 
stead of being relative to the local background stars. 

5.5. Comparison with the USNO-Bl catalog 

Figure 16 compares LSPM proper motions with the proper 
motions quoted in the USNO-Bl.O. More than 75% of the 
stars fall within 30 mas yr"' of the fiusNO-Bi.o = IJ-lspm line- 
The remaining objects are scattered around, with no obvious 
correlation with ^lspm- In particular, there are 975 stars that 
have a USNO-B 1 .0 counte rpart with a p roper motion of ex- 
actly zero. As discussed inl Gould & Kollmeier (2004 ) these 
are not stars with measured proper motions of 0.0, but rather 
stars for which no significant proper motion was calculated (to 
within errors) in the construction of the USNO-Bl.O catalog. 
In any case, these zero-proper-motion stars were identified as 
low-proper motion objects in the USNO-Bl.O. 

We emphasize that the high proper motion status of all 



Fig. 13. — Mean and dispersion in the difference between SUPERBLINK 
and UCAC2 proper motions, as a function of proper motion. The continu- 
ous line shows the mean of the difference, and the dotted lines show the Icr 
dispersion. The accuracy of the SUPERBLINK proper motions is largely in- 
dependent of the magnitude of the proper motions, as demonstrated here. The 
fluctuations above 350 mas yr"' are due to small number statistics. 



LSPM stars has been systematically confirmed by visual in- 
spection. Ambiguous matches with USNO-B 1 .0 counterpart 
have also all been resolved visually. Furthermore, the com- 
parison between LSPM and NLTT proper motions shows a 
very good agreement. Thus, the only explanation for the out- 
liers in Figure 16 is that their USNO-B 1 .0 proper motion is in 
error 

The completeness of the USNO-Bl.O for high proper mo- 
tion stars was investigated by Gould ( 2003b), from a compar- 
ison with the rNLTT. No USNO-Bl.O counterpart was found 
for « 10% of the rNLTT stars ("missing" objects). About 
1 -2% of the matched objects were also found to have dis- 
crepant ("bad") USNO-Bl.O proper motions. The fraction of 
bad or missing stars (dubbed "problem fraction") was found 
to be a function of both magnitude and galactic latitude. The 
reason for this is that it is both more difficult to pair-up detec- 
tions of fast moving stars from different epochs, and to obtain 
accurate proper motions of bright stars, which are saturated 
on the POSS plates. 

We calculate our own "problem fraction" for the USNO- 
Bl.O, using all single stars in the LSPM catalog as a refer- 
ence. We exclude close doubles from the analysis, as they are 
sometimes not resolved in the USNO-B 1 .0, and so might tend 
to overestimate "problem fraction". We calculate the "prob- 
lem fraction" by finding all LSPM stars with no USNO-Bl.O 
counterpart ("missing"), or with a USNO-Bl.O counterpart 
that has WflispM- flusNO-B\.o\\ >40mas yr"' ("bad"). We find 
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Fig. 14. — Difference between the 2MASS-derived 2000.0 positions of 9.105 LSPM stars with the positions of their UCAC2 catalog counterparts. The 
distribution has a dispersion in [RA,DEC] of [91,88] mas. This provides an estimate of the accuracy of the 2000.0 epoch positions of LSPM stars with no 
TYCHO-2 counterparts. 



that most of the problem stars are not "missing" objects, like 
in thelGould (2003b) analysis, but are mostly (> 95%) "bad" 
counterparts. Why this difference? Gould (2003b) matched 
rNLTT stars to USNO-Bl.O objects within a 5" radius, while 
our own search radius was much larger (up to 1')- We have 
thus recovered most of those "missing" USNO-Bl.O counter- 
parts. It appears that those stars have large position errors 
in the USNO-B 1 .0 because they also have large proper mo- 
tion errors. The two values (position and proper motion) are 
linked, because the extrapolated 2000.0 positions are depen- 
dent on a good estimate of the proper motion. In other words, 
the counterparts were "missing" because of their very "bad" 
recorded proper motion. This is why our USNO-B 1 .0 "prob- 
lem" stars are comprised of mostly "bad" counterparts. 

In Figure 17, we plot our USNO-Bl.O problem fraction 
as a function of magnitude, proper motion, and galactic lat- 
itude. A comparison with photographic V magnitude shows 
that large proper motions errors are more common for USNO- 



Bl.O counterparts with 11 < V < 13 and V > 19. The excel- 
lent agreement between the USNO-Bl.O and LSPM for the 
very brightest stars (V < 10) reflects the fact that both the 
LSPM and USNO-Bl.O use TYCHO-2 positions and proper 

motions in t hat rang e. 

Following iGouldl ( l2003b t). we interpret the larger problem 
fraction of the brighter stars from the fact that these objects 
are saturated on the POSS plates, making proper motion de- 
terminations prone to large errors. The large problem fraction 
of faint (y > 19) stars is simply explained by the fact that it 
is much more difficult to pair-up faint objects detected in dif- 
ferent epochs. This increase in the problem fracti on near the 
faint star limit doesn't show up in the Gould ( 2003b) analysis, 
because the analysis was restricted to stars with counterparts 
in the USNO-A catalog, which has a brighter faint magnitude 
limit. 

Large proper motion errors on USNO-Bl.O counterparts 
are also more frequent for stars with larger proper motions 
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Fig. 15. — Comparison between NLTT and LSPM proper motions, for 
all stars in common between the two catalogs. Note that errors on Ulspm 
are < 8" yr"' (see Figure 1 1); the measurement eiTors on ^inltt are cleai'ly 
larger This accounts for the fact that 3 , 000 NLTT stars have LSPM proper 
motions below the fiducial limit of the NLTT catalog (p, = 0.18" yr"'). 



(jJispM > 0.3" yr"'), confirming again the analysis of lGouldl 
(|2003b). It is simply more difficult to pair- up stars that have 
moved very substantially between different photographic sur- 
vey epochs. Finally, crowding is also a major source of confu- 
sion, leading to erroneous USNO-Bl.O proper motions. This 
is evidenced by the increase in the problem fraction at low 
galactic latitude. This is much more significant than the low- 
galactic latitude problem fraction of .Gould c2Q03b), which 
was high only because the low-galactic latitude rNLTT stars 
contain a larger proportion of bright, saturated objects. The 
high problem fraction at low galactic latitude observed by 
iGould (2003b) was thus more the result of those stars be- 
ing dominated by brighter, saturated objects. Our problem 
fraction is a more correct assessment of the completeness of 
the USNO-Bl.O at low galactic latitude: the problem fraction 
reaches 20% near the galactic plane. 

The conclusion is that the USNO-B 1 .0 catalog is at best ap- 
proximately 90% complete and accurate at high galactic lati- 
tude, for stars with fj, < 0.5" yr~' and 14 < V < 19. Otherwise 
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Fig. 16. — Comparison between LSPM and USNO-Bl.O proper motions, 
for LSPM stars that have USNO-Bl.O counterparts. About 75% of the stars 
have USNO-Bl.O proper motions within 0.02" yr"' of LSPM proper mo- 
tions. Several thousand stars have fJ-usNO-Bi.o containing large en'ors. The 
proportion of these"problem" stars, with proper motion errors exceeding 40 
mas yr"' (circle in bottom plot), is significant (see next figure). There is a 
mean offset in the proper motion differences because LSPM lists absolute 
proper motion while USNO-Bl.O lists relative proper motions. 



the completeness falls to 70%. The LSPM is significantly 
more complete in high proper motion stars (see §7 below). 

5.6. SUPERBLINK proper motion accuracy and the 
astrometric magnitude equations 

The SUPERBLINK proper motions are derived from mea- 
surements of photographic Schmidt plates (POSS-I, POSS- 
II). One important issue with these plates is that the mea- 
sured positions of stars have a weak dependence on magni- 
tude. The non-linear response of photographic plates com- 
bined with asymmetric stellar images (from, e.g., imperfec- 
tions in the optical system, or inaccurate guiding) causes the 
photographic image centroids of the stars to be slightly off- 
set compared to that on a linear detector. Because they arise 
from the non-linear response of the substrate, these offsets are 
a function of the magnitude of the star, and they are gener- 
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Fig. 17. — Fraction of stars in the LSPM catalog with no recorded USNO- 
B 1 .0 counterpart, or with a large error in the USNO-B 1 .0 proper motion (as 
compared with the SUPERBLINK proper motion). Close binary stars are 
excluded from the analysis, as they may not be resolved in USNO-B 1.0. The 
calculated USNO-B 1.0 "problem fraction" is lowest for moderately faint stars 
(14 <V< 19) with low proper motions (p, < 0.3" yr"') at high galactic 
latitudes (\h\ > 20), for which it is !a 10%. 



ally largest for bright, saturated objects. These offsets, which 
are also generally dependent on the position on the plate, are 
referred to as the "astrometric magnitude equations". Accu- 
rate astrometric measurements with photographic plates re- 
quire the determi nation a nd application of a magnitude equa- 
tion correction (Girard et al. 1998). 

For the POSS plates, which we use in our SUPERBLINK 
survey, Monet et al. (2003) have shown the existence of fixed- 
pattern astrometric offsets which are on the order of 0. 1 -0.5" 
and are strongly correlated with the XY position on the plate. 
Plotted for stars of different magnitudes, these patterns show 
striking differences, which indicate the existence of signifi- 
cant astrometric magnitude equations. The offsets tend to be 



larger at bright magnitudes, and are also larger near the edges 
of the plates. However, the patterns are relatively regular, and 
much of the variability occurs on scales equivalent to at least 
several minutes of arc. 

As described in §2.1, the SUPERBLINK software intrin- 
sically corrects for large-scale plate distortions by measuring 
the relative proper motions of stars in a frame defined by the 
local background of objects, typically all stars within about 
7' of the high proper motion target. Because SUPERBLINK 
uses the local backdrop of stars to calculate relative proper 
motions, and not the XY plate positions, distortions on the 
photographic plate on scales > 7' do not introduce any signif- 
icant error on the SUPERBLINK proper motions. This how- 
ever, is generally true only if all the stars (target and back- 
ground) are locally offset by the same amount. For example, 
if the local SUPERBLINK reference frame is defined by 16th 
magnitude stars, then the measured relative proper motions of 
16th magnitude stars will not be affected by the astrometric 
magnitude equations (since the target and the reference stars 
are all offset by the same value). This, of course, is true even 
if the value of the offset is different on the first and second 
epoch images. 

On the other hand, if a proper motion target is significantly 
brighter or fainter than the background (reference) stars, the 
astrometric magnitude equations will generally introduce sys- 
tematic errors in the SUPERBLINK proper motions. One ex- 
ception to this case is if the differences in the offsets between 
the target and the reference stars are the same in both the first 
and second epoch images. This may happen, e.g., if the object 
was recorded at the same plate position (X,Y) at both of the 
first and second epoch. Unfortunately, this is generally not the 
case for POSS-I and POSS-II, since their plate centers are on 
different grids. 

Because the systematic proper motion errors introduced by 
the astrometric magnitude equations depend on the difference 
between the magnitude of the target and the magnitude of the 
local background stars, their effect is very difficult to model. 
One would need to estimate the local mean offset of the back- 
ground stars, which have a variable range of magnitudes, and 
use estimates of the offsets expected for the target star given 
its magnitude. One would need to make these estimates sepa- 
rately for the first and second epoch images. Such a procedure 
would be extremely complex. We have thus made no attempt 
to correct the SUPERBLINK proper motions for the effects 
of the astrometric magnitude equations. It is very possible 
that the larger errors on the SUPERBLINK proper motions 
(« 8 mas yr"') compared to the proper motion errors from the 
rNLTT (« 5.5 mas yr~') are due to the fact that astrometric 
magnitude equations have been neglected in SUPERBLINK 
proper motion calculations. 

6. THE CATALOG 
6.1. Format 

The complete catalog in ASCII format is available with the 
electronic edition of this paper The catalog contains 61,977 
lines, each 286 characters long. Each catalog entry consists 
of 29 fields; these are described in Table 1. 

The first field gives the LSPM catalog name. The next 9 
fields provide identifications in the LHS, NLTT, Hipparcos, 
Tycho-2, ASCC-2.5, UCAC-2, 2MASS, USNO-B 1 catalogs, 
when these exist for the star An additional field gives the 
original name of the star in the published literature (e.g. the 
LSR stars of Lepine, Shara, & Rich). In the current version 
of the LSPM catalog, however, the original name is provided 
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TABLE 1 

The LSPM Catalog - field description 



field 


datum 


units 


format 


I 


1 ^ L.? 1 iVJ. UaLulVJ^ IlaillC 




al6 


2 


LHS catalog ID 




a6 


3 


NLTT catalog ID 




a6 


4 


Hipparcos catalog ID 




a7 


5 


Tycho-2 catalog ID 




al2 


6 


ASCC-2.5 catalog ID 




a8 


7 
8 


UCAC-2 catalog ID 
Other name 




a9 
a31 


9 


2MASS catalog ID 




al7 


10 


USNO-B 1 catalog ID 




al3 


11 


R.A. 


degrees 


fl2.6 


12 


Decl. 


degrees 


fll.6 


13 


total relative proper motion 


" yr-' 


f8.3 


14 


relative proper motion in R.A. 


" yr-' 


f8.3 


15 


relative proper motion in Decl. 


" yr-' 


f8.3 


16 


total absolute proper motion 


" yr-' 


f8.3 


17 


absolute proper motion in R.A. 


" vr-' 


f8.3 


18 


absolute proper motion in Decl. 


" yr-' 


f8.3 


19 


astrometric source flag 




a2 


20 


optical B magnitude 


mag 


f6.2 


21 


optical V magnitude 


mag 


f6.2 


22 


photographic blue (Bj) 


mag 


fS.l 


23 


photographic red (Rp) 


mag 


f5.1 


24 


photographic near-IR (In) 


mag 


f5.1 


25 


infra-red J 


mag 


f6.2 


26 


infra-red H 


mag 


f6.2 


27 


infra-red K, 


mag 


f6.2 


28 


estimated V magnitude 


mag 


f7.2 


29 


estimated V—J color 


mag 


f6.2 



only if the star doesn 't have a counterpart in any of the cat- 
alogs listed above. At this point, it is provided as a means 
to distinguish "rediscovered" LSPM stars from those that are 
genuine, "new" discoveries. 

The astrometric information (position, proper motion) is 
detailed in the next 9 fields, and includes a flag that gives the 
origin of the astrometry. The next 8 fields provide the photo- 
metric information, with optical, photographic, and infrared 
magnitudes. The last two fields give the estimated V magni- 
tude and V -J color index. 

A sample of the catalog is shown in Tables 2 and 3, in which 
the first 12 lines are displayed as an example. The full catalog 
is available only in electronic format. 

6.2. Names and identifications 

We assign a LSPM name to each star in our catalog, which 
is based on the star's R.A. and Decl. at 2000.0 epoch in the 
International Celestial Reference System (ICRS, essentially 
equivalent to J2000). The first four characters ("LSPM") are 
the catalog identifiers, and stand for Lepine & Shara Proper 
Motion. A space then separates the catalog ID from the posi- 
tional description. A "J" follows, which indicates the equinox 
of the position. The next four digits are the hours and minutes 
of the R.A., then comes the sign of the Decl. ("H-" for all our 
stars, since we do not have southern declinations - yet) fol- 
lowed by 4 more digits that represent the degrees and minutes 
of Decl. Finally, there is one last character used to distinguish 
stars that would otherwise have the same name. Pairs of stars 
with the same hours/minutes in R.A. and degrees/minutes in 
Decl. have their names appended with an "N", "S", "E", or 
"W" suffix. The choice of suffix depend on the orientation 
of the pair. Their separation in both R.A. and Decl. is deter- 



mined. If the separation in Decl. is larger, then the stars are 
given a N/S suffix, with N ("North") assigned to the star at 
higher declination, and S ("South") to the other star. If it is 
the separation in R.A. that is the largest, then E/W suffixes are 
used, with E ("East") assigned to the star at larger R.A., and 
W ("West") to the other one. By no means are stars with "NS" 
or "EW" suffixes necessarily common proper motion doubles. 
While this is often the case, there are a number of chance 
alignments for which two unrelated high proper motion stars 
happen to be in the same arcminute position bin. Conversely, 
not all common proper motion doubles have "NSEW" suf- 
fixes, since it is often the case that long-period doubles have 
angular separations large enough to put them in separate ar- 
cminute bins. 

Identifications are given for the 2,572 LSPM stars also 
listed in the LHS catalog and for the 31,361 stars listed the 
NLTT. The identification number for stars listed in the LHS 
catalog is their LHS #, which has been traditionally used in 
the literature. The identifications for stars listed in the NLTT 
catalog are the record # in the original NLTT table ("recno" 
in the electronic version of the NLTT catalog, at the VizieR 
catalog service^). 

We also provide identifiers for the 4,839 stars listed in the 
Hipparcos catalog (HIPP number), as well as for 7,943 stars 
with a TYCHO catalog number. A total of 4,306 of the HIPP 
stars also have data in the TYCHO-2 catalog. We give ASCC- 
2.5 identification for a total of 11,430 stars; these include all 
the TYCHO-2 stars as well as the HIPP stars that are not in the 
TYCHO-2. We give UCAC2 catalog numbers for the 9,137 
LSPM stars that are listed in the UCAC2 catalog. Note that 
since the bulk of the LSPM stars are fainter than the magni- 
tude limits of these catalogs, the majority of LSPM stars do 
not have HIPP TYCHO-2, ASCC-2.5 or UCAC2 identifiers. 

A few hundred LSPM stars are not listed in any of the cat- 
alogs listed above, but are not entirely "new" objects because 
they have been previously reported in the literature. An ad- 
ditional column provides the Simbad^ designation for those 
objects. About a third of these are the "LSR" high proper 
motion stars found by our team ( Lepine, Shara, & Rich 2QQ% 
2003). The other stars are a mix of objects, some from old 
proper motion motion catalogs (but that had not been included 
in the NLTT or LHS), some identified as M dwarfs or white 
dwarfs in field spectroscopic surveys. In particular, stars with 
2MASS designations are objects identified recently in various 
searches of ultra-cool M and L dwarfs. 

Finally, we give identifications for the counterparts of the 
LSPM stars in the 2MASS All-Sky Point Source Catalog and 
in the USNO-B 1.0. Catalog identifiers provided in the LSPM 
make it easy to retrieve all relevant information from those 
catalogs. 

6.3. Positions and proper motions 

Positions in the LSPM catalog are given at the 2000.0 
epoch in the ICRS system, and for the great majority are ob- 
tained either by extrapolating from the TYCHO-2 position 
using the TYCHO-2 proper motion (for stars with TYCHO- 
2 counterparts), or by extrapolating from the position of the 
2MASS counterpart, using the SUPERBLINK-derived abso- 
lute proper motion or the ASCC-2.5 proper motion (see below 
for which proper motion is used). 

For stars with no 2MASS counterpart (2,345 objects) we 

* http://vizier.u-strasbg.fr/ 
^ http://simbad.u-strasbg.fr/ 
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used the coordinates calculated by SUPERBLINK from the 
DSS scans instead. The positions of those stars, extrapo- 
lated from the position of the star on the POSS-II scan, are 
much less accurate than the 2MASS -derived positions. Un- 
fortunately, the accuracy of the SUPERBLINK-derived J2000 
positions is only w 0.5" (as estimated from a comparison with 
the 2MASS catalog), and future efforts will be devoted to ob- 
taining more accurate coordinates for those objects. 

In the catalog, we list the relative and absolute proper mo- 
tions in separate columns. The relative proper motions are 
always those determined by SUPERBLINK. This makes it 
easy to identify stars that have not been measured with SU- 
PERBLINK: values for their relative proper motion are set to 
zero. 

For the absolute proper motion, we use either the value 
derived from SUPERBLINK, or quote the absolute proper 
motion from the TYCHO-2 or ASCC-2.5 catalogs. The or- 
der of priority is as follows: (1) TYCHO-2 proper motions, 
(2) SUPERBLINK proper motion, (3) ASCC-2.5 proper mo- 
tion. The quoted TYCHO-2 proper motion errors are < 8 mas 
yr"', for all LSPM stars with a TYCHO-2 counterpart. This 
is smaller than the estimated SUPERBLINK proper motion 
errors, and justifies that we always defer to the TYCHO-2 
proper motion. The quoted proper motion errors for ASCC- 
2.5 stars that do not have TYCHO-2 counterparts are gen- 
erally > 12 mas yr"' with a mean value ~ 14.5 mas yr"', 
larger than the SUPERBLINK eiTors, so we use the ASCC- 
2.5 proper motions only for those stars that have no SU- 
PERBLINK proper motions. The source of the proper mo- 
tion is indicated by the astrometric flag: "T" if the TYCHO-2 
proper motion is used, "S" if it is the SUPERBLINK-derived 
proper motion, and "A" if the ASCC-2.5 proper motion is 
quoted. 

A total of 508 stars in the LSPM catalog have no TYCHO-2 
counterparts nor ASCC-2.5 counterparts, and their proper mo- 
tions have not been measured by SUPERBLINK either For 
those objects, we obtain their proper motion from a variety 
of sources, mainly from the NLTT and rNLTT catalogs, and 
al so from the catalog of revised p roper motions of LHS stars 
of lBak os. Sahu. & Nemeth H2002h . These objects are denoted 
by the astrometric flag "O". 

Note that for LSPM stars with TYCHO-2 counterparts that 
have also been measured with SUPERBLINK, values for the 
absolute proper motion are quoted from the TYCHO-2 cata- 
log, but values for the relative proper motions are those from 
SUPERBLINK. Note also that zonal corrections are calcu- 
lated individually for each star, using all TYCHO-2 objects 
within 7° of that star. Our map of the zonal corrections can 
thus be recovered from the LSPM catalog, by differencing the 
relative and absolute proper motions for stars that have their 
absolute proper motions derived from the SUPERBLINK val- 
ues. 

6.4. Magnitudes 

The LSPM catalog lists optical B, V magnitude from the 
TYCHO-2 and ASCC-2.5 catalogs. It also gives the photo- 
graphic blue (Bj), red (Rp), and near infrared (I^) magnitudes 
extracted from the USNO-B 1 .0 catalog. Stars with no USNO- 
Bl counterparts are listed with Bj and In magnitudes of 99.9, 
indicating these to be unavailable. When a value for Rp could 
not be obtained from the USNO-B LO catalog, we used the 
value estimated by SUPERBLINK from the POSS-II red DSS 
scans. A value of 99.9 is also used whenever there is only par- 
tial magnitude information from the USNO-B 1 counterpart 



(or counterparts, see §3.8 above). 

Infrared magnitudes for LSPM stars are obtained from 
their counterparts in the 2MASS All-Sky Point Source Cat- 
alog (Cutri et al. 2003). The accuracy is about 0.02 mag for 
5 < 7 < 14, 5 < // < 14, 4 < < 13; it is w 0.25 mags for 
brighter stars (saturated on the 2MASS images). For fainter 
stars, the uncertainty increases with magnitude. The 2MASS 
catalog is complete to J~16.5, H~16.0, and Kj ~ 15.5 ^. 

Finally, estimated V magnitudes and V -J colors are given 
in the last two columns. Values are provided for all but a few 
entries. These should be used for quick reference, and for 
classification of the high proper motion stars. See §4.4 for the 
caveats in using these estimated values. 

7. COMPLETENESS OF THE LSPM CATALOG 

7.1. Comparison with the NLTT catalog 

The completeness of a proper motion catalog can be a func- 
tion of magnitude, proper motion, and position. The com- 
pleteness is generally dependent on how easy it is to identify 
moving objects against the backdrop of the "fixed" stars. De- 
tection will be more difficult if the star is fainter, but also if it 
is moving faster, or if the local density of background objects 
is larger Traditionally, proper motion surveys have been very 
incomplete for faint stars at low Galactic latitudes. As de- 
scribed above, our SUPERBLINK software, with its image 
subtraction algorithm, was specifically designed to address 
this problem, and detect moving stars in densely populated 
areas. We therefore expect the LSPM catalog to be signifi- 
cantly more complete than the NLTT at low galactic latitudes. 
The main question is how much more complete the LSPM is. 
In particular, we would like to know whether the LSPM cat- 
alog still suffers from some incompleteness near the Galactic 
plane. 

We first compare the distribution of NLTT and LSPM stars 
as a function of position, separating the stars into two groups: 
stars brighter than V = 16, and stars fainter than V = 16. Fig- 
ure 18 shows the distribution of brighter stars. The distribu- 
tions for both catalogs are relatively uniform (no sharp dis- 
continuities, or "holes"). However, the density of objects still 
appears to be non-uniform, and is larger in high Galactic lati- 
tudes. However, the decline is very gradual, in sharp contrast 
to the distributions of faint NLTT stars, shown in Figure 19. 
The density of faint NLTT objects falls very abruptly at low 
Galactic latitudes: a clear mark of true incompleteness in the 
NLTT catalog. On the other hand, the distribution of faint 
LSPM objects follows more or less the same trends as the dis- 
tribution of bright high proper motion stars. 

Two interpretations are possible. One is to say that the dis- 
tribution of high proper motion stars should be uniform across 
the sky, and that it is the completeness of the catalog that pro- 
gressively diminishes as one goes to lower galactic latitudes. 
With this interpretation, both the NLTT and LSPM are signif- 
icantly incomplete, even at moderately high Galactic latitudes 
(20 < l^?! < 60). The other interpretation, and the one we fa- 
vor, is that the distribution of proper motion selected objects 
is naturally non-uniform over the sky, and that the progressive 
density variations observed in Figure 18, and in the top panel 
of Figure 19 have little to do with completeness. Under this 
second interpretation, the completeness of the LSPM catalog 
is uniformly high, both at high and low galactic latitude. We 
demonstrate the truth of this statement as follows. 

Detailed documentation can be found at 
littp://www.ipac. caltech.edu/2mass/ 
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Fig. 18. — Distribution of high proper motion stars from our LSPM catalog (top) compared to the distribution of stars from the NLTT (bottom). Shown here 
are stars brighter than V = 16.0 (see Figure 19 for the distribution of fainter stars). The Galactic equator is shown (thick line). Note that even though both catalogs 
are expected to be relatively complete (> 90%) in that range, it is obvious that the density of objects is larger at high galactic latitudes. This suggests that proper 
motion selected samples are intrinsically non-uniform. Indeed, one does expect proper motion surveys to be more sensitive to old disk and halo stars at high 
galactic, because of tiie asymmetric drift. 



Figure 20 lets us estimate the completeness of the LSPM 
catalog from the rate at which NLTT stars are recovered by 
SUPERBLINK and TYCHO-2. One can guess in advance 
that the rate is very high, since very few NLTT stars had to 
be separately incorporated into the catalog (see §3.4). We ex- 
clude from this analysis the NLTT stars that were missed by 
SUPERBLINK because they are in areas thai were not pro- 
cessed by the code, since their inclusion would underestimate 
the true efficiency of SUPERBLINK. We calculate the frac- 
tion of NLTT stars that have been recovered either by SU- 
PERBLINK or from the TYCHO-2/ASCC-2.5 catalogs. Re- 
sults are shown in Figure 20. The recovery rate is « 99% to 
a magnitude as faint as V = 19, falling to w 90% at y = 20. 
Note also the small dip (to « 98%) at the boundary between 



TYCHO-2 and SUPERBLINK stars (V = 12), which we in- 
vestigate further in §7.2. The recovery does not vary signifi- 
cantly with the proper motion. There is however the expected 
trend that SUPERBLINK misses more stars at low galactic 
latitudes, but the recovery rate still exceeds 97% at \b\ < 10. 

The recovery rate of NLTT stars by 
SUPERBLINK/TYCHO-2 is a good estimate of the 
completeness of the LSPM catalog, but is valid only in 
regions of the [iJ,V,b] parameter space that contain a suffi- 
cient number of NLTT stars. We know from Figure 19 that 
the NLTT catalog is significantly depleted at \b\ < 15 and 
V > 16, which means the completeness of the SUPERBLINK 
sample cannot be evaluated for faint stars at low Galactic 
latitudes using the above method. Outside of that specific 
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Fig. 19. — Same as Figure 18, this time showing stars fainter than V = 16.0. The NLTT is dramatically incomplete in low-galactic latitudes. This situation is 
much improved with the LSPM catalog. One still observes a lower density of high proper motion stars at low galactic latitudes. 



range, however, we conclude that the completeness of the 
SUPERBLINK/TYCHO-2 sample is indeed extremely high. 

But the completeness of the LSPM catalog itself is larger 
than the completeness of the SUPERBLINK/TYCHO-2 sam- 
ple, since the missing NLTT stars have been included in the 
LSPM. Assuming that the NLTT catalog is itself more than 
90% complete for > 15 and V < 18 stars suggests that 90% 
of the stars missed by TYCHO-2 and SUPERBLINK might 
have been found by Luyten, in which case the LSPM catalog 
could up to w 99.9% complete down to V = 18. However, 
this assumes that both the SUPERBLINK/TYCHO-2 and the 
NLTT samples are statistically independent, an assumption 
that may not be entirely valid. Indeed, it is very possible that 
stars that have been missed by SUPERBLINK have also been 
missed by Luyten for the exact same reason. From our expe- 
rience, failed detections mostly occur when a faint star moves 
on pixels saturated by a bright nearby objects at one of the 
two epochs (or worse, at both epochs). But high proper mo- 



tion stars eventually move out of the glare, so a star that is 
hidden at one epoch will be easy to spot at another Indeed, 
this is what happens for several of the NLTT stars that SU- 
PERBLINK failed to identify: the star is in the glare of a 
bright object on the POSS-II image. Because Luyten used his 
own 1960s plates as a second epoch , the object wa s then easy 
to identify. As was shown in Le pine. Shara. & Rich (,2003i) . 
many of the new high proper motion stars found with SU- 
PERBLINK follow the inverse pattern: easy to find on the 
POSS-II scans, their trajectory puts then in the glare of a 
brighter star in the 1960s. The real problem is for stars hid- 
ing on the POSS-I plates, since both we and Luyten are using 
it as the first epoch. Additionally, there might be a few faint 
stars lost in the extended, saturated patches from very bright 
stars at all epochs (if the total 45 years motion of the star is 
shorter that the size of the saturated image), although these 
cases should be quite rare. 
The bottom line is that for every star missed by SU- 
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Fig. 20. — Estimated completeness of the LSPM catalog based on the 
recovery rate of NLTT stars. The recovery rate is > 98% for all values of the 
proper motion and for magnitudes V > 19.0. The recovery rate then drops to 
~ 90% at y = 20.0. The recovery rate is marginally smaller at low galactic 
latitudes, but still exceeds 97%. 



Fig . 21. — Completeness of the old NLTT catalog, estimated from the frac- 
tion of LSPM catalog objects that are NLTT stars. The NLTT completeness 
is a function of both magnitude and proper motion. While the NLTT was 
relatively complete for stars with fi > 0.3" yr"' and brighter than Ri.=l4, 
its completeness dropped significantly at fainter magnitudes, and was much 
more limited for proper motions near the cutoff at ;i = 0.18" yr~^ . 



PERBLINK but found by Luyten, there is probably at least 
another one that has been missed by both. Thus, despite the 
fact that the additional NLTT stars make the LSPM catalog 
more complete than the SUPERBLINK/rYCHO-2 sample, 
we conclude that the LSPM catalog is approximately 99.0% 
complete to as faint as V = 19, andathigh (\b\ > 15) Galactic 
latitudes. 

We now determine the completeness of the NLTT catalog 
by comparing it to the LSPM catalog, assumed to be 99% 
complete within the limits quoted above. We first plot the 
fraction of LSPM stars that are NLTT objects as a function of 
magnitude, proper motion, and galactic latitude (Figure 21). 
The general completeness of the NLTT over the whole north- 



ern sky for stars with fi > 0.15" is a little above 60% down to 
y = 14, falls gradually to 40% at V = 19, and then drops more 
abruptly. This is, however, not a fair assessment of the com- 
pleteness of Luyten's survey, because the LSPM catalog has 
a lower proper motion limit, and because the completeness of 
the NLTT is significantly lower at low galactic latitudes. 

First we want to determine the effective proper motion limit 
of the NLTT, specifically the limit above which the NLTT is 
most complete. We find that the completeness of the NLTT 
catalog decreases sharply below /i = 0.20" yr"'. This is to be 
expected from fiducial limit of 0.18" yr"^ and the 0.02" yr"' 
measurement error. Because the drop in completeness starts 
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Fig. 22. — Completeness as a function of magnitude for the NLTT catalog, 
based on a comparison with the LSPM, for various ranges of proper motion 
and galactic latitudes. The "Completed Palomar Region" is traditionally as- 
sumed to be the most complete (top), although the NLTT is significantly more 
complete if we exclude the bands 10 < \b\ < 15 (middle). At low Galactic 
latitudes (bottom) the completeness of the NLTT drops significantly beyond 
y = 15. 



a little above fi = 0.20" yr~', we define the range of maximal 
completeness of the NLTT as /i > 0.25" yr"'. 

The NLTT is supposed to be most complete within the 
Completed Palomar Region (CPR), as defined by Dawson 
{1986), which for the northern sky is simply \b\ > 10, for 
stars with fi > 0.2" yr"'. However, we note that the complete- 
ness of the NLTT actually starts falling a few degrees above 
\b\ = 10. We conclude that the NLTT is most complete for 
proper motions /i > 0.25" yr"' and galactic latitudes \b\ > 15. 

We now proceed to check the completeness of the NLTT as 
a function of magnitude for three regions of parameter space 
(Figure 22). First we check the completeness for the CPR. We 



find that the estimate of lDawsoiil Jl986h is accurate, and that 
the NLTT catalog is indeed more than 80% complete down 
to y = 18. The second region is the one we identified as the 
most complete of the NLTT, with the ^ > 0.25" yr"' stars at 
Galactic latitudes \b\ > 15. For that restricted area, we find 
the NLTT to be 90% complete down to V = 18.5. Finally, we 
check the completeness of the NLTT at low Galactic latitudes 
1^1 < 15, again for stars with /i > 0.25" yr"'. We now find 
that the completeness falls from about 90% at V = 15.0 to 
only 30%belowy = 17.0. 

Our estimate of a fairly high completeness of the NLTT at 
high Galactic latitudes is very significant, because is indicates 
that the internal completeness test described bv IFlvnn et al\ 
(I^OOl) underestimates the completeness of the proper mo- 
tion sample. The test of'F lvnn et al\ (2QQX) suggested that the 
completeness of the NLTT in the CPR fell to 80% at V = 15 
and down to 60% at V = 18.5, which is significantly smaller 
than the results of our (external) test. It appears that the crit- 
icism offered by Monet ef a/. (2000) is legitimate, and that 
changes in the space density of objects as a function of dis- 
tance do lead to an undere stimate of the comp leteness when 
applying the internal test of lFlvnn 'erall(l200!l) . 

7.2. Completeness at the SUPERBLINK/TYCHO-2 boundary 

Completeness problems in the LSPM catalog occur in 
the magnitude overlap region between TYCHO-2 and SU- 
PERBLINK stars (V ~ 12). The problem arises because the 
completeness of the TYCHO-2 catalog starts to decrease be- 
fore SUPERBLINK reaches its full detection efficiency. As a 
result, it is very possible that relatively bright stars have been 
missed by SUPERBLINK because of plate saturation, and at 
the same time are missing from the TYCHO-2 (or ASCC-2.5) 
catalogs because they are fainter than these catalogs' com- 
pleteness limits. 

Fortunately, the TYCHO-2 catalog has a small, but signifi- 
cant overlap with the SUPERBLINK detections. A large frac- 
tion of the brighter SUPERBLINK stars are in the TYCHO-2 
catalog and, likewise, a large fraction of the fainter TYCHO- 
2 stars have been recovered by SUPERBLINK. There is 
no significant magnitude gap in the combined TYCHO-2 
-I- SUPERBLINK sample. This is clear from Figure 20 
(top panel), which demonstrates that the combined TYCHO- 
2/SUPERBLINK sample does recover more than 98% of the 
NLTT stars in that magnitude range. Since the completeness 
of the full LSPM (including the missed NLTT stars) is very 
high, we can use the assumption that it is in fact complete 
to estimate the completeness of both the TYCHO-2 and SU- 
PERBLINK samples as a function of magnitude. The esti- 
mated SUPERBLINK and TYCHO-2 completeness is plotted 
in Figure 23 (top two panels). 

We use a probabilistic approach to estimate the complete- 
ness of the combined sample. We can say that the com- 
pleteness function f{V) is also the probability that a star of 
V magnitude will be in the sample. The probability that a 
star is in the TYCHO-2 (or ASCC-2.5) catalog is friV) ~ 
Nt(V)/Nlspm(V), while the probability for a star to be de- 
tected by SUPERBLINK is fs(V) ^ Ns(V)/Nlspm(V), where 
N(V) is the number of star in the sample as a function of 
magnitude. The probability that a star is missing from the 
TYCHO-2 is thus 1 -fr, while the probability that a star will 
be missed by SUPERBLINK is 1 -fs- Assuming the two sam- 
ples to be independent, we can say that the probability that a 
star is missing from both the TYCHO-2 and SUPERBLINK 
samples is (l-frX^-fs)- Thus, the completeness of the com- 
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Fig. 23. — Estimated completeness of the LSPM catalog at the magnitude 
boundary between the TYCHO-2/ASCC-2.5 catalog and the SUPERBLINK 
stars. The LSPM is buUt from the combination of the two samples, which 
overlap over a relatively small magnitude range. The coverage appear to 
be tight enough that few stars should have been missed at the boundary; the 
estimated completeness of the combined sample only drops to 98% around 
V = 12.0, compared to si 99.5% at brighter magnitudes, and 99% at fainter 
magnitudes. The estimated number of stars missing from the LSPM because 
of the boundary effect should be < 100. Note that several of these were 
recovered as additional NLTT objects (see Figure 7). 



bined sample will be 1 - (1 -/r)(l -/s), which is a function 
of V magnitude. From this, we calculate the completeness 
of the combined SUPEBLINK+TYCHO-2 list, and the total 
number of stars that would presumably be missing from the 
LSPM because of this incompleteness (bottom two panels in 
Figure 23). The effect is very small, and the completeness 
only drops to « 98% around V = 12. Overall, we expect to be 
missing w 50-100 stars because of the boundary effect. 

We note that « 60 additional NLTT stars have been found in 
thel0<y<14 magnitude range that were neither in the SU- 
PERBLINK nor TYCHO-2/ASCC-2.5 hsts (Figure 7). The 



NLTT thus recovers some of the missing stars. But because 
the NLTT has a higher proper motion limit than the LSPM, 
additional stars in the 0.15" yr"' < ^ < 0.18" yr"' range are 
probably still missing. Since about a third of the LSPM stars 
are in the 0.15" yr~^ < /* < 0.18" yr"' range, we are probably 
still missing « 30 stars in the 10 < V < 14 magnitude range. 

7.3. Completeness of the LSPM at low galactic latitudes 

A close examination of Figures 18-19 makes it very appar- 
ent that the number density of LSPM stars is not uniform over 
the sky. Clearly, there are more LSPM stars at higher galac- 
tic latitude than near the Galactic plane. The lower counts at 
low galactic latitude are actually a natural consequence of the 
low proper motion cut-off of the catalog, and not a result of 
decreased completeness in low-galactic latitude fields. 

The first line of evidence that this is true is that the lower 
density of LSPM stars at low galactic latitude is observed 
both for bright and fainter stars. Figure 25 shows that there 
are ~ 40% fewer LSPM stars at low Galactic latitudes than 
there are at high Galactic latitudes, and this is independent of 
the magnitude of the stars. If there were completeness prob- 
lems in the LSPM because of crowding, or other low galactic 
latitude effects, the proportion of low Galactic latitude stars 
would drop with fainter magnitudes, as is observed in the 
NLTT catalog (see Figure 24). 

If the distribution of LSPM stars is not uniform over the 
northern sky, it must be because of proper motion selection 
effects. The velocity components of the stars in the vicinity 
of the Sun are not isotropic, because of Galactic rotation and 
because of the Sun's motion relative to the local standard of 
rest. The distribution of stellar proper motion vectors is very 
much dependent on the position on the sky, as illustrated in 
Figure 25. One naturally expects to find more high proper 
motion stars at high Galactic latitude because of the large ap- 
parent drift of the halo and old-disk stars in that direction, as 
seen from the Sun. 

If the distribution of high proper motion stars were uniform 
over the sky, estimating the completeness at low galactic lat- 
itudes would be straightforward. Since the LSPM is nearly 
complete at high Galactic latitudes (for V < 19), the complete- 
ness at low Galactic latitude could then be estimated from the 
ratio between the number density of low Galactic latitude stars 
to the number density of high Galactic latitude stars. This 
could be calculated for each magnitude bin to obtain the com- 
pleteness as a function of magnitude. However, because the 
distribution of high proper motion stars is not uniform over 
the sky, one has to predict what the density of stars at low 
Galactic latitude should be. Estimates of the completeness 
will be dependent on the predicted number density of high 
proper motion stars in the low Galactic latitude regions. 

To estimate the completeness of the LSPM in low galactic 
latitude fields, we use the following test. For stars in some 
magnitude bin [V, V + AV], we first calculate the number den- 
sity of LSPM objects (in stars per square degrees) in the area 
located within 15 degrees of the Galactic plane {\b\ < 15). 
We then calculate the number density of [V, V + AV] stars in 
10 degree bands located above and below this region (15 < 
\b\ < 25). We assume that the survey is complete in the 
15 < \b\ < 25 region, and that the distribution of high proper 
motion stars should be uniform over the whole \b\ < 25 area. 
The completeness in the \b\ < 15 region is thus simply the ra- 
tio between the measured density in \b\ < 15 to the measured 
density in 15 < |^| < 25. We repeat the calculation for a range 
of magnitude bins, to obtain the completeness as a function of 
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Fig. 24. — Surface density of (in stars per square degrees) of tlie NLTT and LSPM catalogs as a function of galactic latitude for three groups of /i > 0.2" yr 
stars: bright (top), moderately faint (middle) and faint (bottom). At high galactic latitudes, outside of the area delimited by the dashed line (|fe| > 20), the NLTT 
is > 90% complete, and the LSPM > 99% complete, for V < 18. Both catalogs are > 90% complete for bright stars (V < 14) in low galactic latitudes {\b\ < 20). 
The drop in the density of V < 14 NLTT stars at low Galactic latitudes is largely due to catalog incompleteness. Not so for the LSPM catalog, which appears to 
be largely complete at < 20, given the intrinsic trend in number density with b. 



V. Note that this method essentially provides an internal test 
of completeness. The main caveat is that the density of high 
proper motion stars is not uniform over the sky, as demon- 
strated above. However, the use of relatively low Galactic 
latitude bands (15 < \b\ < 25) as a reference should mini- 
mize the effects of the intrinsic non-uniformity. Nevertheless, 
since the density of high proper motion stars decreases with 
Galactioc latitudes, our internal completeness test is expected 
to sUghtly underestimate the actual completeness level at low 
Galactic latitude, possibly by as much as « 5 - 10%, although 
a detailed modeling of the distribution of high proper motion 
stars would be required to determine the exact value. 

Results of the completeness test are shown in Figure 26, 
where we have apphed it separately to the complete LSPM 
catalog (top), and for the subsample of LSPM stars that are 
in the NLTT catalog (bottom). Results for the NLTT stars are 



consistent with the external completeness test shown in Fig- 
ure 22 (which is based on a comparison between NLTT and 
LSPM), with a sharp drop in completeness at V = 16. Note 
that at moderately bright magnitudes (13 < V < 15) the in- 
ternal test suggests that the NLTT is only « 80-85% com- 
plete, which is lower than estimated from the external test 
(« 95%, see Figure 22). There are two different interpreta- 
tions for the differences in the two NLTT completeness es- 
timates. First, if the LSPM is only 90-95% complete at 
< 15, 13 < y < 15 (as suggested by the internal test), then 
the external test overestimates the NLTT completeness, be- 
cause it assumes the LSPM to be 100% complete. If indeed 
the LSPM is only 90-95% complete, then the external test 
possibly overestimates the completeness by 5 - 10%, which 
would bring both values in closer aggrement. The second in- 
terpretation is that the external test is right, and the LSPM is 
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Fig. 25. — Di,stribution of absolute (zonal corrected) proper motions for LSPM stars within 15° of the north Galactic pole (top), the direction of Galactic 
rotation (bottom right), and the Galactic anti-center (bottom left). The distribution of proper motion vectors is not uniform, and shows a distinct pattern that is 
clearly dependent on the position on the sky. The fixed, low-/i cutoff of the catalog (empty disks centered on the origin) determines how many stars are locally 
included in the LSPM catalog. The larger scatter and larger offset in the proper motion vectors near the Galactic pole implies that more stars in that direction will 
make it into the LSPM. This explains the higher density of LSPM objects at high galactic latitudes, near the north galactic pole. 



100% at \b\ < 15, 13 < y < 15, but the internal test underes- 
timates the completeness level because there are intrinsically 
fewer stars at \b\ < 15 than at 15 < \b\ < 15. If the second 
interpretation is valid, then the internal test indeed underesti- 
mates the completeness by as much as 10%. 

The internal test suggests that the LSPM catalog is at least 
80% complete for stars brighter than V = 19. This should be 
taken as a the lowest possible value; as discussed above, the 
internal test may be underestimating the completeness by up 
to 10%. Thus, it is possible that the LSPM is actually « 90% 
in low Galactic latitude regions. For stars fainter than V = 19, 
the internal completeness test shows a drop to « 40% just be- 
low V = 20. The LSPM catalog is thus definitely shallower at 
low Galactic latitudes, by 1 to 2 magnitudes. While the LSPM 
should be regarded as largely complete to at least V = 20.0 
at high Galactic latitudes, it should be regarded as relatively 
complete only down to V = 19.0 in low Galactic latitude re- 



gions. 

Our completeness estimates at low Galactic latitudes are 
relatively crude at this point, and we are sorry we cannot pro- 
vide more accurate values. More accurate completeness esti- 
mates could be obtained only if we had better estimates of the 
expected density of high proper motion stars in the low Galac- 
tic latitude regions. Alternatively, a proper modeling of all 
the effects that result in SUPERBLINK missing high proper 
motion stars might also provide more accurate completeness 
estimates. We believe the first option to be beyond the scope 
of this paper The second option, in our opinion, would be ex- 
tremely difficult to carry out, because of the complexity of the 
SUPERBLINK algorithm, but also because one would need 
to characterize the combined effects of all the pairs of POSS- 
I/POSS-II plates used, with their variety of saturation levels, 
point spread functions, and a complete assessment of plate 
defects. Ultimately, the most accurate estimates of the com- 
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Fig. 26. — Completeness at low galactic latitude, estimated from the ratio 
between the number density of high proper motion stars at \b\ < 15 and the 
number density of high proper motion stars at 15 < \b\ < 25. This internal 
test suggests that the LSPM is > 80% complete at low Galactic latitudes for 
V < 19 (top). Applied to NLTT stars (bottom) this completeness test yields 
results that are very similar with the completeness estimate shown in Figure 
22, with a sharp drop at V = 16. This internal test probably underestimates 
the completeness by 5 — 15%, since there are intrinsically fewer high proper 
motion stars at lower Galactic latitudes (see Figure 24). 



pleteness levels will come from better, more accurate proper 
motion surveys at low Galactic latitudes (even of limited ar- 
eas), which will provide true external tests to the complete- 
ness of the LSPM catalog. 

8. STELLAR CONTENTS OF THE LSPM CATALOG 

8.1. Color-magnitude classification 

The apparent V magnitude of LSPM stars is plotted in Fig- 
ure 27 as a function of the V -J color index. The LSPM stars 
are nicely clumped in three main groups. The interpretation is 
straightforward when one considers the selection effects im- 



plied in a sample of stars selected for high proper motions. 
Catalogs of high proper motions stars essentially contains 
subsamples of nearby stars. However, the spatial extent of 
the detection range depends on the transverse velocity of the 
star. Hence, high velocity (halo, old-disk) stars tend to be se- 
lected from a larger distance. It is thus fair to say that, on first 
approximation, high proper motion catalogs combine nearby 
disk stars with more distant halo/old-disk stars. Because of 
the limiting distance, a diagram of apparent magnitude as a 
function of color will have the same general features as a 
color-magnitude diagram: a main sequence extending from 
upper left to lower right, and a white dwarf sequence at the 
bottom left. This is indeed what we see in Figure 27. Because 
the stars are not exactly at the same distance, all sequences are 
thicker and fuzzier, but they are recognizable. Because high 
velocity stars are detected to a larger distance, they will form 
their own sequence, but shifted down because or their larger 
distance and hence larger mean apparent magnitude. Again 
this feature is quite apparent in Figure 27. 

Note in Figure 27 the concentration of stars on the lower 
left of the plot, where the nearby white dwarfs are expected 
to be found. All stars with a 2MASS counterpart are plotted 
in red, while stars with no 2MASS counterparts are plotted in 
blue. It is clear that the majority of the high proper motion 
stars with no 2MASS counterparts are white dwarfs, which 
confirms the conjecture posed by Salim & Gould ( 2003). 

In the H/J-K color magnitude diagram, the brighter stars 
are separated into two main groups: a dense clump of stars 
around J-K5=0.8, and a more diffuse one extending from the 
blue edge of the first one to around J-K5=0.3 (Figure 28). 
While the blue clump extends only down to about H=13, the 
red clump extends all the way down to the magnitude limit, 
where the larger magnitude errors scatter the stars about in 
J-Ks- The big red clump is populated by M dwarfs and sub- 
dwarfs, which are degenerate in J-K color from M2 to M7 
(|Bessell & Brett 1988). The blue, diffuse clump consists of 
F-G-K dwarfs and subdwarfs. The fact that the distribution of 
main sequence stars other than M dwarfs ends (at H=13) well 
before the 2MASS magnitude limit indicates that the LSPM is 
complete for those stellar subtypes. Note also the diffuse wisp 
that extends from [J-K5,H]=[0.9,7.5] to [L5,4.5]; these are the 
very few giant stars that have proper motions large enough to 
be included in the LSPM. 

Dwarfs and subdwarfs are expected to occupy distinct loci 
on the color magnitude diagram because these populations 
have different mean distances in proper motion catalogs. The 
proper motion limit (/i > 0.15" yr"') restricts the detection 
range of disk objects, whose transverse velocities are typi- 
cally < 50 km s"', while halo objects (subdwarfs) with typ- 
ical transverse velocities 2-5 times as large, can find their 
way into the proper motion catalog from distances 2-5 times 
larger This is illustrated in Figure 28, where we have noted 
the probable positions of subdwarfs, 3-4 magnitudes below 
the dwarfs. 

We also note in Figure 28 the locations expected to be pop- 
ulated by white dwarfs, and late-type M dwarfs (M7-M9) and 
L dwarfs. Of course, one has to look for those objects above 
the limit (H=15) below which 2MASS colors become less re- 
liable. We nevertheless find a significant number of object that 
are convincing white dwarf candidates. Late-type M dwarfs 
(M7-M9) ai-e found in the range 1 .0 <J- Kt < 1.3, while L 
dwarfs normally occur beyond J-K, > 1.3 ("Kirkp atrick et al\ 
1999) . Several candidates are detected, and they warrant 
further investigation. While many might be known objects. 
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Fig. 27. — Optical/infrai'ed color-magnitude diagram of all stars in the LSPM catalog. Stars with TYCHO-2 of ASCC-2.5 counteipails are plotted in green; 
these have the most accurate values of V and V-J. Stars with no 2MASS counterparts are plotted in blue; most of them are consistent with being white dwarfs. 



chances are high that at least a few will be new ones, in partic- 
ular if they occur at low galactic latitudes, which have largely 
been avoided by previo us L dwarf surveys (iKirkpatrick et al\ 
fT999tlCruz et aJ2003l) . 

The J-H/H-K color-color diagram is largely consistent with 
the vast majority of LSPM objects being main sequence stars 
with no significant reddening (Figure 29). The distribution 
is very similar to the general color-color diagram of 2MASS 
sources in high galactic latitude fields, and follows the stan- 
dard sequence determined by Bessell & Brett ( 1988). The loci 
of different types of objects are indicated, as in Figure 28. 

Overall, the accuracy of the 2MASS J, H, and Kj mag- 
nitudes allows one to separate different object subtypes and 
populations (much better than the photographic magnitudes). 
They should be very useful in planning follow-up observa- 
tions of the LSPM stars. 



8.2. Reduced proper motions 



Reduced proper motion diagrams are a major tool in the 
classification of loc al stars into different stellar populations 
tToneslll972t lEvansi ri992: Salim & Gould 2000). It was re- 
cently demonstrated by S alim & Gould (|2002) that optical- 
infrared reduced proper motion diagrams are very efficient in 
separating samples of high proper motion stars into three dis- 
tinct classes: main sequence disk dwarfs, halo subdwarfs, and 
white dwarfs. 

Using the proper motions and magnitudes in the LSPM cat- 
alog, we build a reduced proper motion diagram for all but 27 
LSPM stars using a reduced proper motion calculated from 
the V band {Hy) and the V-J color. The reduced proper 
motion is analogous to an absolute magnitude in which the 
proper motion is used in place of the parallax. While the ab- 
solute magnitude is defined as 

My =y + 5*log7r + 5, 

where tt is the parallax in seconds of arc, the proper motion is 
defined as 

i/y = y + 5*log/l + 5, 
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0.5 1 

Fig. 28. — Infrared color-magnitude diagram of stars in our proper motion catalog, based on JHKj magnitudes extracted from the 2MASS All-Sky Point 
Source Catalog. Main sequence F-G-K stars populate the upper left of the distribution, while M dwarfs are clumped around J-Kj = 0.8. A diffuse wisp of giant 
stars can be seen at the upper right. Below H=15, errors in 2MASS magnitudes increase, which blurs out the color distribution. 



where is the proper motion in seconds of arc per year. The 
reduced proper motion is directly related to the absolute mag- 
nitude: 

Hv=Mv + 5 *\ogVT -3.38 

where vr is the projected velocity of the star in the plane of 
the sky, in km s~' . 

The reduced proper motion diagram is thus analogous to 
a color-magnitude diagram, except for the fact that the usual 
stellar sequences are "blurred" by the logvj- term. But dif- 
ferent populations have different ranges of possible vj. Disk 
stars, on the one hand, have a mean transverse velocity (vr) ^ 
50km s~\ which yields: 

{Hv)disk = Mv + 5.l. 

Halo stars, on the other hand, have significantly larger mean 
transverse velocities (vj-) ~ 300km s"\ so their Hr are gener- 
ally larger: 

{Hv)halo=Mv + 9.0. 



The (disk) dwarfs and (halo) subdwarfs thus form distinct se- 
quences in the reduced proper motion diagram, with halo sub- 
dwarfs located well below the disk dwarfs. White dwarfs oc- 
cupy a distinct locus on the lower left of the diagram, a posi- 
tion familiar to users of color-magnitude diagrams. 

We show in Figure 30 the reduced proper motion diagram 
for LSPM stars. The corrected, absolute proper motions were 
used. The loci of the different stellar classes and popula- 
tions are labeled. Note the similarity with Figure 27, which 
plotted the apparent magnitude as a function of color. The 
two plots are fundamentally very different, however. The re- 
duced proper motion is a function of luminosity and velocity, 
while the apparent magnitude is a function of luminosity and 
distance. The two plots look similar only because stars are 
piled up against the LSPM proper- motion limit of 0. 15" yr~\ 
which introduces a correlation between the velocity and dis- 
tance. High-velocity stars can make it into the catalog even 
at large distances, while only the nearest of the low-velocity 
stars have proper motions large enough to be in the LSPM. 
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Fig. 29. — Infrared J-H/H-Kj color-color diagram of stars in the proper motion catalog. The distribution is consistent with stars on the main sequence with no 
reddening. 



Figure 27 thus separates the halo stars based on their larger 
average distances, while Figure 30 separates stars based on 
their larger transverse velocities. For classification purposes, 
the reduced proper motion diagram should be preferred. 

The reduced proper motion diagram contains a wealth of 
information, and the separation between the different popula- 
tions warrants a more detailed analysis of the individual pop- 
ulations represented here. Such a detailed analysis is beyond 
the scope of this paper. Nevertheless, we wish to emphasizes 
the significant potential of the LSPM catalog in the study of 
the local stellar populations. 

9. CONCLUSIONS 

We have generated a new catalog of stars with proper mo- 
tions larger than 0.15" yr"' which currently is the most com- 
plete of its kind. We have achieved our initial goals of locat- 
ing new high proper motion stars, and redetermining to higher 
accuracy the positions and proper motion estimates of previ- 
ously known objects, especially the high proper motion stars 
from the LHS and NLTT catalogs. 

The catalog is Umited to northern declinations, and lists 



61,977 stars down to a magnitude V = 21. This essentially 
doubles the number of previously cataloged high proper mo- 
tion stars in this hemisphere. The catalog is estimated to be 
> 99% for y < 19 stars at high Galactic latitude {\b\ > 15), 
and « 90% for V < 19 stars at low Galactic latitude (\b\ < 15). 
This is a very significant improvement over previous catalogs. 

We provide photometric estimates in the three optical bands 
of the POSS-n survey: Bj (HlaJ emulsion with GG385 filter), 
Rf (IIIaF emulsion with RG610 filter), and % (IVN emul- 
sion with RG9 filter). We also provide B and V magnitudes 
for the brighter stars, from their TYCHO-2/ASCC-2.5 cata- 
log counterparts. While the B and V magnitudes are very ac- 
curate, much room for improvement remains with the optical 
photometry of fainter stars. Infrared photometry, on the other 
hand, was extracted from the 2MASS All-Sky Point Source 
catalog, and is accurate down to about 15th magnitude in each 
band (J, H, Kj). The probability of our stars having been mis- 
matched in the 2MASS catalog is almost zero. 

To compare all stars in the same color/magnitude system, 
we provide an estimate of the apparent V magnitude and the 
V -J color index for all but 814 LSPM stars. These are most 
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reliable for the brighter (V < 12) stars, but should be used 
with caution for fainter stars, for which V is estimated from 
the USNO-Bl.O photographic magnitudes. 

The LSPM catalog is a work in progress. An extension to 
the southern sky is currently underway. We are also working 
on an expansion to lower proper motions, down to 0. 10" yr~' . 
Future plans include an improvement of the magnitudes, es- 
pecially in the optical bands, by using magnitude estimates 
from a variety of other sources (SDSS, USN0-A2 catalog, 
GSC-2.2 catalogs, future versions of the UCAC catalog). 

At this point, the LSPM catalog is ideally suited for follow- 
up observations of selected targets of interest. Indeed, we are 
currently working on a massive spectroscopic survey of se- 
lected LSPM objects, including all stars with proper motions 
/i > 0.45" yr"' , for which spectroscopic observations are now 
in hand (Lepine et al. 2005, in preparation). 

The LSPM catalog will be updated as new discoveries are 
made. We invite investigators who discover new high proper 
motion objects in the northern sky to contact the authors so 
that their discovery can be included in the catalog. Likewise, 



investigators who notice that a known high proper motion star 
is missing from the LSPM, or who find errors in our data are 
invited to conmiunicate with us. 
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APPENDIX 

NLTT STARS NOT INCLUDED IN THE LSPM-NORTH CATALOG 

A total of 1,859 objects listed in the NLTT catalog, and claimed to be located north of the celestial equator, did not make it into 
the LSPM catalog. The reasons why these stars were not included in the LSPM catalog fall into four categories: (1) the object is 
moving slower than the 0.15" yr"' Umit of the LSPM catalog, (2) the object does not show up on the POSS plates, (3) the object 
is a duplicate NLTT entry, (4) the object is a high proper motion nebula. 

Most of the slow moving objects are stars that are incorrectly listed in the NLTT as having large proper motions. A total of 208 
such objects are stars that are listed in the TYCHO-2 catalog; their TYCHO-2 proper motions unambiguously place them under 
the LSPM inclusion Umit. An additional 11 04 stars had their proper motions remeasured with SUPERB LINK, and the updated 
value makes them low proper motion (/x <0.15" yr~') stars. 

There are also 76 objects which have a quoted NLTT proper motion under the 0.15" yr~^ limit of the LSPM. While some of 
these are possibly bogus (see below), 4 of them have their low proper motion status confirmed in the TYCHO-2 catalog. Another 
63 have their low proper motion confirmed by SUPERBLINK. Three more stars did not have their proper motion remeasured by 
SUPERBLINK but direct examination of DSS scans confirms they are low proper motion stars. 

We have identified 39 objects to be dupUcate entries of other NLTT stars. In all cases, the dupUcate entry was Usted with a 
slightly different position, generally within a few arcminutes of the primary entry. 

The remaining 395 objects could not be recovered on the POSS plates, and are listed as "bogus". An NLTT star can be missing 
if it initially was a spurious detection. This is highly probable for stars initially reported to have a magnitude near the plate limit 
(V>18). Another possibility is that the quoted NLTT catalog position was very far off from the actual location of the star. In this 
case, it is very likely that the object has been picked up by SUPERBLINK and is now listed in the LSPM catalog as a "new" high 
proper motion star. Thus the object is not really missing, but was rather lost, and has now been rediscovered. 

Two high proper motion "stars" in the NLTT are found to be small, compact nebulae with large proper motions. The two 
objects are NLTT 13414 and NLTT 13424. Moving nebula are not included in the LSPM catalog at this point. These, and other 
candidate high proper motion nebulae found with SUPERBLINK will be discussed in an upcoming paper. 

All NLTT objects north of the celestial equator that are not in the LSPM catalog are listed in the accompanying table (see 
Table A-1 for a sample of the first 15 lines). The table lists the NLTT catalog number, followed by the position, red and blue 
magnitude, and proper motion as quoted in the NLTT catalog. The NLTT positions are transformed into J2000 coordinates from 
their original values. A flag value of "B" indicates that the star could not be found by direct examination of the Digitized Sky 
Survey scans of the Palomar Sky Survey (POSS) photographic plates. A flag value of "D" means that the object is a duplicate 
NLTT catalog entry. A flag value of "L" means that the star has an NLTT proper motion under the nominal lower limit (0.15" 
per year) of the LSPM catalog. A flag value of "N" means that although the object was found to be real, it is not a star but rather 
a proper motion nebula. A flag value of "S" means that the proper motion of the star was remeasured by SUPERBLINK, and 
found to be under the nominal lower Umit (0.15" per year) of the LSPM catalog. A flag value of "T" means that the star is in the 
TYCHO-2 catalog, and that its TYCHO-2 proper motion is under the nominal lower limit (0.15" per year) of the LSPM catalog. 
The corrected position and proper motion is given for all stars that could be recovered. 

We note that this list does not contain a single star from the LHS catalog. Every one of the northern LHS stars has been 
accounted for and is now included in the LSPM catalog. 



